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ABSTRACT 

We present a detailed study of the largest sample of intervening O vi systems in the red- 
shift range 1.9 < z < 3.1 detected in high resolution (R ~ 45,000) spectra of 18 bright QSOs 
observed with VLT/UVES. Based on Voigt profile and apparent optical depth analysis we find 
that (i) the Doppler parameters of the O vi absorption are usually broader than those of C iv 
(ii) the column density distribution of O vi is steeper than that of C iv (iii) line spread {6v) 
of the Ovi and Civ are strongly correlated (at 5.3cr level) with 6v(Ovi) being systematically 
larger than Sv{C iv) and (iv) 6v(0 vi) and 6v(C iv) are also correlated (at > 5cr level) with their 
respective column densities and with A^(Hi) (3 and 4.5cr respectively). The median column 
densities of H i, O vi, and C iv are found to be higher when low ions are present. A^(C iv) and 
A^(H i) are strongly correlated (at 4.3cr level). However, no significant correlation is found be- 
tween A^(0 vi) and A^(H i). These findings favor the idea that C iv and O vi absorption originate 
from different phases of a correlated structure and systems with large velocity spread are prob- 
ably associated with overdense regions. The velocity offset between optical depth weighted 
redshifts of C iv and O vi absorption is found to be in the range < |Av(0 vi - C iv)| < 48 
km s"' with a median value of 8 km s"'. 

We do not find any evidence for the ratios A^(0 vi)/A^(H i), N(Ovi)/N(Crv) and 
A^(C iv)/A^(H i) to evolve with z over the redshift range considered here. But a lack of sys- 
tems with high A^(0 vi)/A^(H i) ratio (i.e., > -0.5 dex) for z > 2.5 is noticeable. Similar trend 
is also seen for the A^(Civ)/A^(Hi) ratio. We compare the properties of Ovi systems in our 
sample with that of low redshift (z < 0.5) samples from the literature and find that (i) the O vi 
components at low-z are systematically wider than at high-z with an enhanced non-thermal 
contribution to their /^-parameter, (ii) the slope of the column density distribution functions for 
high and low-z are consistent, (iii) the range in gas temperature estimated from a subsample of 
well aligned absorbers is similar at both high and low-z, and (iv) Qq yj = ( 1 .0 + 0.2) x 1 0"^ for 
A^(0 vi) > lO'-'' cm"^, estimated in our high-z sample, is very similar to low-z estimations. 
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1 INTRODUCTION 

The study of low density intergalactic medium (IGM) is ex- 
tremely important because it forms the primary reservoir of baryons 
throughout the cosmic ages. These baryons get accumulated into 
galaxies in the process of structure formation. The heavy elements 
produced in galaxies got transported to the IGM by means of out- 
flows driven by supemovae or tidal interactions. Thus the IGM en- 
richment history provides useful constraints on the star formation 
history and contribution of various feedback mechanisms at dif- 
ferent epochs. The tenuous IGM is detectable in the form of Lyo- 
and heavy element absorption lines in the QSO spectra. Hence the 
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observation of Lyo- and metal lines are crucial to understand the 
interaction between galaxies and the surrounding IGM. 

The observations of Civ and Ovi absorption in QSO spec- 
tra have establis hed the presence of he avy elements in th e IGM 
unequivocally (Cowieetal. 'l995a"'^; Songaila & Cowie 'l99d; 
Bergeron et al. 2002; Simcoe et al. 2002; Carswell et al. 2002|). 
The early observations of heavy elements in the high redshift uni- 
verse were focused on C iv absorption, since it has strongest lines 
falling in the region free from Lya contamination. Several high red- 
shift surveys have ascertained that the cosmic density of C iv ab- 
sorb ers has no t evolved subs tantially from redshift z = 5 to z = 1.5 
(see ISongailj [20 01. 2005; iBoksenberg et all |2003| ; IPettini et all 
l2003l ; ISchave et al..20 03). Due to lack of sufficient QSO sight lines, 
a clear picture is yet to emerge regarding the evolution of C iv at 
z > 6. However the limited sample of Civ absorbers available 
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until now indeed indicates that a good fraction of metal s may al- 
ready be present in the IGM even at these high redshifts I Songailal 
200d : ISimcoj|2006l : iRvan-Weber et alJlioM : ISecker et alJbOOd 



20091) . The enrichment level found by these studies is consistent 
with [C/H] ~ -2.8. 

Given the low metallicity, the direct detection of metals in the 
underdense regions [with overdensity, 6{= nn/nH) ^ 10], which 
occupy most of the volume of the universe at any given epoch, 



methods like pixel analysis are used instead jEUison et alj 


200C 


Schaye et alj 2003 


: Aracil et al."2004': 'Aauirre et alj 2001 


2008 


Scannapieco et alj 


2006; Fieri et al. 2006). They show that metals 



must be present even in the underdense regions. However the frac- 
tional volume occupied by the metals is still unknown. 

Even in regions where metal absorption is detected directly, it 
is unclear what are the main physical processes that maintain the 
ionization state of the gas. In general, it is believed that photoion- 
ization by the meta-galactic UV background keeps the gas ionized. 
On the other hand, the winds that seed the IGM with metals and 
the accretion shocks in the evolving density fields may also pro- 
vide sufficient mechanical feedback to collisionally ionize the gas. 
Therefore, it is crucial to simultaneously study different species 
covering a wide range of ionization states to get a better understand- 
ing of the metal enrichment and the different ionizing mechanisms 
at play. 

Under photoionization by UV background, the Ovi absorp- 
tion is generally produced from regions of low-density having high 
ionization parameter. In addition, the high cosmic abundance of 
oxygen makes O vi a good tracer of metal-enrichment in the low- 
density IGM. In fact photoionization seems to b e a viable pro- 
cess f or most of the high r e dshift O vi absorbers jBergeron et al.l 



I2OO2I : ICarswell et d] I2OO2I : [Bergeron & Herbert-Fort 200! ). On 



the other hand, hy 


Dave et al.' 


1999, 


Smith et al. 


2011; 



'2001'; 'Fang & Bryan' '200 1'; iKang et alj 



T999I ; 



20051 ; 



amount of baryons could reside in the warm-hot phase of the inter- 
galactic medium (called WHIM with T x 10^ - 10' K) and this frac- 
tion evolves with redshift. Highly ionized species of oxygen such as 
O VI, O VII and O viii can be useful probes of the WHIM. In the op- 
tical regime, O vi is the best species to probe relatively cooler phase 
(i.e., r ~ 3 X 10^ K) of the WHIM because, the ionization fraction 
of O VI has its maximum around this temperature in case of colli- 
sional ionization. Indeed, it has been suggested that a large fraction 
of th e O VI absorption associated with the IGM ( Sim coe et al .I2OO2I 
I2OO6) and high-z damped Lya systems (F ox et al, 2007h may origi- 
nate from collisionally ionized gas. Thus the origin of the ionization 
of the intergalactic O vi absorbers is a matter of debate. 

Numerous extragalactic Ovi systems have been detected 
at z < 0.5 with the Far Ultraviolet Spectroscopi c Explorer 
(FUSE) and Hubble Space Telescope {HST) (e.g., 
2000; Stockeetal. 2006; Danforth et al. 2006; iTric 



Tripp et akl 
etalj|2008l; 



Thom & Chen 2008 ^Wakker & Savageii2009l ; iLehner et alj|2009l ; 
Tumlinson et al. 201 ij). However these studies could not establish 
convincingly whether O vi absorbers are predominantly tracing the 
WHIM gas or not. They seem to arise in eithe r high metallicity 



photoionized (e.g., lOppenheimer & Dave 



ity collisionally ionized (e.g., ISmith et al 



20091) or low metallic 
20 III) gas loosely as- 



sociated with galaxies. An alternative tool to detect the WHIM, 
which is independent of the metal enrich ment, is to search for 
thermally-broadened Lya absorbers (BLA s) ('Semb ach et al.|[2004l ; 
Richter et 'Zll2004l2006l ;l Danforth et al ][2010; Sava ge et 



also be en ambiguous. The few detections of Nevm in the FUV 
regime jSavage et alJl2005l ; lNaravanan et al.ll2009ll20Ilh and O vn 



20Ilh . However, until now the WHIM detection through BLAs has 



in the soft X-ray regime (e.g..'Nicastroetal."2005a|lbl-] Buote et all 
[20 09 ; Fang et al. 2010; Zappacosta et al. 2010) possibly hint the 
existence of the WHIM. 

At high re dshift (z > 2 ) previ ou s studies jSimcoe et alj 

| 2002l |2004. ICarswell et all |2002|; [Bergeron et alj |2002|; 

'Bergeron & Herbert-Fort' '2005'; ISchave et alj |2007| ; iFrank et alj 
i2010a b) have already provided important insights into the prop- 
erties of Ovi absorbers. Here we present a detailed analysis of 
Ovi absorbers using a sample which is twice as large as the 
previous sample of high-z intervening O vi systems studied at high 
resolution. 

This paper is organized as follows. In section[2]we describe the 
observations and the data reduction procedure for our data sample. 
In section[3]we describe the line identification strategy and present 
the O VI and C iv sample and various physically motivated subsam- 
ples. In section|4]we analyze the distributions of Ovi Voigt profile 
parameters and compared them with those of C iv. In this section 
we also compare the properties of O vi absorption at high and low 
redshift. In section [5] we discuss the line kinematics of Ovi and 
C IV absorption. In section|6]we present analysis based on total col- 
umn densities with results of photoionization model as guidelines. 
In section|7]we summarize our results. 

Throughout this paper we use the following cosmological pa- 
rameters for a flat universe : = 0.3, Q.\ = 0.7, Cli^h- = 0.02 and 
Ho = 71 km s^'Mpc"'. The solar relative abundances are taken to 
be default values used in CLOUDY v(07.02), i.e., log (C/H)o = 
-3.61 and log (0/H)o = -3.31. 



2 OBSERVATIONS 

The spectra used in this study were obtained with the Ultra- Violet 
and Visible Echelle Spectrograph (UVES) (Dekkeretal. 200^ 
mounted on the ESO Kueyen 8.2 m telescope at the Paranal ob- 
servatory in the course of the ESO-VLT large progr amme "The 
Cosmic Evolution of the IGM" ( iBergeron et alj|2004|) . This large 
programme provided a homogeneous set of 18 QSO sight lines with 
QSO emission redshifts ranging fro m 2.1 to 3.3. The raw data were 
reduced using the UVES pipeline ( iBallester et"ai]|2000l) which is 
available as a dedicated context of the MIDAS data reduction soft- 
ware. The main function of the pipeline is to perform a precise 
inter-order background subtraction for science frames and master 
flat fields, and to apply an optimal extraction to retrieve the object 
signal, rejecting cosmic ray impacts and performing sky subtrac- 
tion at the same time. The reduction is checked step-by-step. Wave- 
lengths are corrected to vacuum-heliocentric valu es using standard 
conversion equations tedledl 1 966l ; Istumpfil 1 980h . Combination of 
individual exposures is performed by adjusting the flux in each in- 
dividual exposures to the same level and inverse variance weighting 
the flux in each pixel. Great care was taken in computing the error 
spectrum while combining the individual exposures. Our final error 
in each pixel is the quadratic sum of the weighted mean of errors 
in the different spectra and the scatter in the individual flux mea- 
surements. Errors in individual pixels obtained by this method are 
consistent with the rms dispersion around the best fitted continuum 
in regions free of absorption lines. The final combined spectrum 
covers the wavelength range of 3000 to 10, 000 A with occasional 
narrow gaps in the red. During the observations, the 2x2 binning 
mode was used yielding a binned pixel size of 2.0 - 2.4 km s"'. 
A typical signal-to-noise ratio (S/N) ~ 30 - 40 and 60 - 70 per 
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Table 1. Details of O vi systems 



QSO 


im 


Zmin 


Zmax 




logWCHi) 


log N{0 vi) 


log N{C IV) 


Class' 


Case- 




iSi<0 VI) 
(kms"') 


(5v(C IV) 
(kms-') 


|Ai'{Ovi-Civ)| 
(kms-') 


Low Ions'" 


HE 1341-1020 


2.135 


1.983 


2.083 


1.9982 


13.77+0.04 


13.71+0.04 


12.22+0.05 


del 


A 





24.3 


19.7 


4.6 


No 










2.0414 


15.70+0.07 


14.01+0.12 


13.25+0.09 


bcl 


B 


+1 


88.6 


36.9 




Yes 










2.0850 


15.10+0.02 


13.79+0.08 


12.98+0.06 


bb 


B 


1 


23.6 


91.8 




No 


Q 0122-380 


2.190 


1.989 


2.137 


2.0349 


15.54+0.06 


13.70+0.09 


13.10+0.04 


del 


B 


1 


59.8 


77.7 




No 










2.0626 


12.49+0.04 


13.47+0.05 


12.94+0.01 


bd 


A 





45.1 


19.3 


2.0 


No 


PKS 1448-232 


2.220 


2.008 


2.166 


2.1099 


13.86+0.02 


14.48+0.07 


13.16+0.03 


dd 


A 





56.0 


35.4 


3.3 


No 










2.1660 


15.52+0.36 


14.29+0.13 


13.48+0.18 


dd 


A/B 





272.8 


166.4 


13.4 


Yes 


PKS 0237-23 


2.222 


1.954 


2.168 


1.9878-* 


13.64+0.01 


13.41+0.06 


< 11.97 


dd 


A 





41.3 


0.0 




No 










2.0108'' 


14.51+0.06 


13.15+0.05 


< 11.93 


dd 


A 





26.5 


0.0 




No 










2.0412 


16.24+0.04 


13.16+0.07 


12.61+0.06 


bd 


B 


1 


31.0 


36.8 




No 










2.0422 


< 13.48 


14.30+0.07 


13.66+0.05 


bd 


A 





47.7 


40.7 


7.2 


No 


HE 0001-2340 


2.263 


1.982 


2.209 


2.0323 


16.07+0.16 


13.89+0.06 


12.85+0.13 


bd 


B 





50.3 


64.2 


15.3 


No 










2.1617 


15.61+0.22 


13.82+0.06 


12.76+0.06 


dd 


B 





80.4 


80.7 


19.2 


No 


Q 0109-3518 


2.404 


1.980 


2.347 


2.0226 


14.99+0.05 


14.24+0.14 


12.52+0.03 


bd 


B 


+1 


60.1 


41.0 




Yes 










2.1415 


15.01+0.53 


13.52+0.06 


11.70+0.08 


dd 


B 





60.1 


16.6 


5.3 


No 


HE 2217-2818 


2.414 


1.959 


2.357 


2.0160 


15.03+0.02 


13.77+0.13 


12.80+0.05 


bb 


B 





86.8 


86.0 


0.1 


No 










2.0748 


14.12+0.01 


14.32+0.01 


12.85+0.01 


dd 


A/B 


-1 


37.8 


92.0 




No 










2.1808 


16.06+0.04 


14.19+0.23 


13.53+0.07 


dd 


A/B 





178.2 


115.1 


46.4 


Yes 


Q 0329-385 


2.435 


2.018 


2.378 


2.0764 


13.70+0.08 


13.26+0.05 


13.21+0.01 


bd 


A 





21.3 


19.2 


4.2 


No 










2.1470* 


14.73+0.21 


13.95+0.18 


< 12.89 


dd 


A 





120.1 


0.0 




No 










2.2489* 


13.51+0.30 


14.25+0.49 


< 12.82 


dd 


A 





107.4 


0.0 




No 










2.2510 


15.82+0.55 


14.86+0.16 


14.41+0.08 


dd 


B 





127.4 


84.4 


10.4 


Yes 










2.3139 


14.23+0.10 


13.30+0.10 


12.35+0.12 


dd 


B 





48.2 


43.1 


12.9 


No 










2.3520 


13.05+0.01 


14.09+0.09 


13.50+0.08 


dd 


A 





41.2 


27.3 


1.0 


Yes 










2.3639 


14.85+0.09 


13.73+0.03 


12.40+0.10 


dd 


A 





64.8 


57.8 


1.2 


No 










2.3738 


15.26+0.01 


14.26+0.08 


12.66+0.06 


dd 


A 





153.0 


96.6 


0.8 


No 


HE 1158-1843 


2.449 


1.980 


2.391 


2.2354 


14.74+0.02 


13.71+0.17 


12.99+0.12 


dd 


A 





47.1 


35.4 


2.1 


No 










2.2660 


15.88+0.05 


13.68+0.05 


13.68+0.04 


bd 


A/B 


+1 


142.4 


29.7 




Yes 


HE 1347-2457 


2.611 


1.985 


2.551 


2.1162 


15.16+0.07 


14.57+0.07 


13.43+0.07 


bd 


A/B 


+1 


91.2 


63.2 




Yes 










2.2349 


15.06+0.02 


14.43+0.47 


13.24+0.11 


dd 


A 





195.9 


187.3 


10.6 


No 










2.3289 


16.45+0.04 


14.63+0.04 


14.26+0.04 


bb 


B 


-1 


59.1 


99.5 




Yes 










2.3327* 


13.89+0.02 


13.26+0.03 


< 12.02 


dd 


A 





45.9 


0.0 




No 










2.3422* 


14.82+0.30 


13.45+0.09 


< 12.30 


dd 


A 





56.7 


0.0 




No 










2.3700 


14.96+0.10 


14.08+0.05 


12.95+0.08 


bd 


A 





77.8 


62.7 


16.9 


Yes 










2.4455 


15.14+0.47 


13.69+0.26 


12.62+0.17 


bd 


B 





67.7 


64.7 


0.8 


No 


Q 0453-423 


2.658 


2.001 


2.597 


2.1694 


13.75+0.03 


14.27+0.05 


12.65+0.05 


bb 


A 


+1 


119.2 


80.5 




No 










2.2765 


16.29+0.06 


14.87+0.32 


14.78+0.23 


bb 


B 


-1 


179.6 


149.8 




Yes 










2.3978 


15.14+0.09 


14.97+0.36 


14.47+0.18 


bb 


A 





186.0 


156.2 


48.5 


Yes 










2.4435 


15.37+0.05 


14.44+0.24 


13.76+0.35 


bb 


A/B 





181.4 


142.8 


6.9 


Yes 










2.5028 


15.77+0.02 


13.99+0.03 


13.71+0.34 


bb 


B 


1 


74.7 


78.3 




Yes 










2.5214 


15.31+0.04 


13.58+0.19 


12.58+0.15 


bb 


B 


-1 


51.6 


32.5 




No 










2.5371* 


14.67+0.08 


13.42+0.07 


< 12.29 


dd 


A 





45.2 


0.0 




No 


PKS 0329-255 


2.703 


2.080 


2.641 


2.2044 


15.44+0.22 


14.61+0.24 


13.32+0.21 


bd 


B 


+1 


176.9 


144.7 




No 










2.3284 


16.20+0.12 


13.96+0.30 


13.49+0.70 


bb 


B 


-1 


93.9 


186.8 




Yes 










2.3896* 


14.22+0.12 


13.37+0.02 


< 12.00 


bb 


A 


+1 


0.0 


0.0 




No 










2.4253 


15.58+0.05 


13.86+0.12 


13.25+0.19 


dd 


B 





65.8 


50.0 


9.7 


Yes 










2.5687 


14.78+0.04 


13.58+0.44 


12.52+0.26 


dd 


A 





42.8 


48.0 


6.8 


No 


Q 0002-422 


2.767 


2.064 


2.704 


2.1767* 


13.62+0.01 


13.72+0.04 


< 12.02 


bd 


A 





52.6 


0.0 




No 










2.2200* 


14.95+0.02 


13.80+0.32 


< 12.09 


dd 


A/B 





51.8 


0.0 




No 










2.2621 


15.47+0.04 


13.67+0.11 


12.43+0.03 


bd 


B 





44.5 


38.5 


12.0 


No 










2.5395 


14.55+0.03 


13.79+0.01 


12.48+0.02 


dd 


A 





61.5 


29.1 


7.2 


No 










2.6075 


14.72+0.03 


14.10+0.05 


12.14+0.05 


bd 


A 





70.4 


34.9 


7.7 


No 


HE 0151-4326 


2.789 


2.043 


2.726 


2.0885 


15.16+0.06 


13.75+0.10 


12.81+0.29 


dd 


B 





35.3 


33.8 


2.7 


No 










2.1699 


15.25+0.02 


13.86+0.07 


13.04+0.01 


bd 


B 





43.5 


27.5 


9.7 


No 










2.2010 


15.39+0.03 


14.23+0.13 


13.12+0.19 


bd 


B 





131.6 


146.9 


1.1 


Yes 










2.3598* 


14.03+0.03 


13.45+0.10 


< 12.24 


bd 


A 





47.6 


0.0 




No 










2.4681 


13.25+0.22 


14.38+0.54 


13.43+0.06 


dd 


A 





73.3 


65.9 


15.0 


Yes 










2.4927 


14.77+0.01 


14.16+0.27 


13.51+0.03 


bd 


A 





120.6 


42.6 


18.1 


No 










2.5053 


15.01+0.02 


14.19+0.08 


12.47+0.16 


bb 


A 


1 


118.1 


117.5 




No 










2.5235 


14.81+0.56 


13.94+0.23 


12.30+0.15 


bb 


A 


-1 


0.0 


0.0 




No 


HE 2347-4342 


2.871 


2.301 


2.806 


2.3475 


16.17+0.11 


13.74+0.11 


13.49+0.01 


bd 


B 





34.7 


25.6 


9.7 


Yes 










2.4382 


14.84+0.02 


14.38+0.11 


12.98+0.15 


bd 


B 





126.8 


74.7 


11.0 


No 










2.6346 


14.84+0.38 


14.31+0.28 


13.08+0.14 


bb 


A 





126.0 


64.4 


26.5 


No 










2.6498 


15.12+0.06 


14.07+0.08 


12.66+0.06 


bd 


A/B 





114.8 


79.8 


2.1 


No 










2.7121 


14.85+0.07 


13.96+0.06 


12.27+0.11 


bd 


B 


-1 


86.8 


46.2 




No 










2.7356 


16.50+0.28 


14.41+0.11 


14.08+0.10' 


bb 


B 


-1 


58.2 


0.0 




Yes 










2.7456 


14.42+0.04 


13.47+0.04 


7 


bd 


A 





71.0 


0.0 




No 


HE 0940-1050 


3.084 


2.458 


3.016 


2.5167 


15.37+0.03 


13.89+0.01 


12.84+0.07 


bb 


B 





62.2 


74.8 


16.6 


No 










2.6433 


15.56+0.15 


13.94+0.07 


13.16+0.07 


bb 


B 


-1 


90.9 


89.5 




No 










2.6580 


15.41+0.02 


13.79+0.01 


13.40+0.05 


dd 


B 





68.5 


32.8 


5.1 


Yes 










2.8265 


16.72+0.37 


14.50+0.06 


14.17+0.16 


bb 


B 


-1 


136.6 


377.0 




Yes 










2.8345 


16.66+0.06 


14.40+0.01 


14.27+0.24 


bb 


B 


-1 


83.2 


91.4 




Yes 










2.9377 


14.77+0.03 


13.86+0.12 


12.94+0.05 


bb 


A 


-1 


101.3 


70.2 




Yes 










2.9401 


14.67+0.01 


14.16+0.06 


12.64+0.09 


bb 


A/B 


+1 


92.2 


71.9 




Yes 


Q 0420-388 


3.117 


2.625 


3.048 


2.6615 


15.44+0.03 


14.58+0.03 


14.37+0.18 


bb 


B 


-1 


229.2 


89.6 




Yes 










2.8100 


15.29+0.03 


13.48+0.05 


12.53+0.04 


bd 


B 





46.3 


17.3 


5.1 


No 


PKS 2126-158 


3.280 


2.508 


3.209 


2.9073 


16.16+0.04 


13.92+0.17 


13.63+0.23 


bd 


B 





68.7 


82.5 


7.6 


Yes 



Table Notes - Underlined lines of siglit are common to Bergeron & Herbert-Fort 1 2005 ). ' Class is based on the O vi profiles. The systems with both the doublets are unblended (partially blended ) are 

marked by 'V^/" {"bb"). In "btf^ systems one of the doublets is blended. ^ Case is based on the presence of unsaturated Lyman-series line. "Case A" systems have at least one of the available 
Lyman-series line unsaturated. In "Case B'\ systems all the available Lyman-series are saturated. "Case A/B" are the cases where some parts of H i absorption are unsaturated. ^ Siypc is based on the 
robust measurement of O vi line spread. 6iypc of + 1 and - 1 indicate upper and lower limits on Hne spread whereas 6iypc of represents robust measurement. "O vi only" system. C iv column 
density is taken from iAeafonova et alj l200'7h . ^ Indicating the presence of low ions. ^ C iv falls in the spectral gap. 
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Table 2. Systems with upper limits on A'(0 vi) 







log (column den.sity) 




5v{C iv) 


Low 


QSO 




^(H I) 


iV(Civ) 


< ^(0 vi) 




Ions 


PKS 1448-232 


1.9516 


14.94 + 0.05 


12.89 + 0.01 


13.90° 


23.9 


No 


PKS 1448-232 


1.9781 


15.07 ± 0.07 


12.64 + 0.03 


13.77° 


35.6 


No 


HE 0001-2340 


2.1634 


14.71 ±0.04 


11.92 + 0.06 


12.81° 


20.2 


No 


Q 0109-3518 


2.0463 


16.09 + 041 


14.07 + 0.10 


14.38 


172.3 


Yes 


HE 2217-2818 


2.0374 


15.46 ± 0.07 


12.16 + 0.06 


12.91° 


23.3 


No 


HE 2217-2818 


2.1553 


14.13+0.02 


12.48 + 0.02 


12.94 


23.9 


No 


HE 1158-1843 


2.0348 


15.43+0.29 


12.61 ± 0.03 


13.28° 


38.8 


No 


HE 1158-1843 


2.0407 


15.52 + 0.07 


12.51+0.03 


13.17° 


29.1 


No 


HE 1347-2457 


1.9750 


14.76 + 0.02 


12.64 ± 0.04 


13.60° 


45.7 


No 


Q 0453^23 


2.4163 


15.00 + 0.01 


12.56 ± 0.02 


13.81 


28.5 


No 


PKS 0329-255 


2.1611 


15.95 + 0.09 


1241 + 0.05 


13.23 


29.4 


Yes 


PKS 0329-255 


2.2953 


14.82 + 0.03 


11.96 + 0.09 


12.57° 


19.1 


No 


PKS 0329-255 


2.4208 


14.86 + 0.02 


12.53 + 0.16 


14.31 


41.7 


Yes 


PKS 0329-255 


2.5868 


15.14 + 0.02 


12.29 ± 0.04 


12.76° 


15.9 


Yes 


Q 0002-422 


2.3647 


12.28 + 0.02 


12.14 + 0.03 


12.56° 


15.3 


No 


HE015M326 


2.4013 


15.13+0.01 


12.56 + 0.08 


13.97 


60.5 


Yes 


HE015M326 


24158 


13.35 + 0.01 


13.03 ± 0.01 


14.19 


23.4 


Yes 


HE 0151^326 


2.4196 


13.04 + 0.02 


12.75 ± 0.01 


12.57° 


21.7 


Yes 


HE015M326 


2.5199 


15.40 ± 0.02 


12.28 ± 0.03 


12.56° 


25.9 


Yes 


HE 2347^342 


2.3132 


15.89 + 071 


13.73 ± 0.07 


14.82 


194.8 


Yes 


HE 2347^342 


2.3317 


15.58 + 0.06 


12.38 ± 0.04 


14.31 


35.9 


Yes 


HE 0940-1050 


2.3307 


16.32 + 0.43 


14.85 ± 0.23 


15.19 


197.2 


Yes 


HE 0940-1050 


2.4090 


15.94 + 0.14 


13.56 + 0.16 


14.18 


35.2 


Yes 


HE 0940-1050 


2.6136 


15.32 + 0.02 


12.56 ± 0.02 


13.58 


31.6 


No 


HE 0940-1050 


2.6679 


15.60 + 0.06 


13.78 + 0.09 


14.11 


68.6 


Yes 


Q 0420-388 


2.8235 


15.52 + 0.03 


13.73 + 0.10 


14.17 


1324 


Yes 


Q 0420-388 


2.8496 


14.25 + 0.04 


12.70 + 0.13 


14.06 


64.9 


No 


Q 0420-388 


2.9519 


15.35+0.01 


12.72 + 0.10 


13.93 


58.2 


Yes 


PKS 2126-158 


2.3889 


13.92 + 0.02 


13.05 + 0.01 


13.66 


19.7 


No 


PKS 2126-158 


2.5537 


14.05 + 0.04 


13.10 + 0.05 


14.55 


28.2 


No 


PKS 2126-158 


2.6790 


14.00 + 0.03 


14.24 ± 0.03 


14.24 


33.9 


Yes 


PKS 2126-158 


2.8194 


15.61 +0.02 


13.50 + 0.08 


14.60 


170.5 


No 


PKS 2126-158 


2.9634 


15.85 + 0.03 


13.27 + 0.09 


14.08 


91.3 


Yes 



Table Note - °Unabsoibed continuum is seen at least in one the O vi doublets. 



pixel was achieved at 3300 and 5500 A respectively. The typical 
final spectral resolution is R ~ 45, 000 (FWHM ~ 6.6 km s"') over 
the entire wavelength range. This spectral resolution allows us to 
resolve lines with fc-parameter as narrow as ~ 4 km s"'. The unab- 
sorbed QSO continuum is then fitted using low order polynomials 
extrapolated from wavelength ranges devoid of strong absorption 
lines. The detailed des cription of data cal ibration is presented in 
lAracil et al.l ( l2004h and lChand et ai] ( l2004h . 



Table 3. List of Lyman limit systems with O vi detections 



QSO 


^sys 


log ^(Hi) 


log ^(0 VI) 


log N{C IV) 


PKS 0329-255 


2.1584 


17.53 


12.99 


12.86 


Q 0002^22 


2.1683 


18.46 


13.64 


14.38 




2.3023 


17.75 


15.21 


15.47 




2.4633 


18.49 


14.63 


14.79 


PKS 2126-158 


2.6377 


19.23 


14.85 


a 




2.7692 


18.62 


14.80 


14.71 



Table Note - " C iv falls in the spectral gap 



Q 0453-423 || z,^, = 2.5371 





C IV 1550 




C IV 154B 

1 . 1 . 1 . 1 . 1 


1 ' 1 ' 
..^■'■u.-^r-^.ii. OJIVV- 




. r VI 1037 

^ ' ' ' ' 








VI 1031 /\ 






^^^^^^ 


H 1 972 




1 ' 1 ' — 1 — 1 ' 1 ' rrz — 


■ yvH 1 1025 /■'^ \ 




: 1 ' 1 '-n— 1 ' 1 ' 1 


H 1 12 15"'"''^^^^^^,^ 




-200 -100 


100 200 



Relative Velocity (km s"') 

Figure 1. An example of O vi system without detectable associated C iv and 
other heavy elements absorption. The best fitted Voigt profiles of O vi and 
H I are over-plotted on the observed data. The vertical tick marks show the 
positions of the individual components. The absoiption redshift that defines 
the zero velocity and the name of the background QSO are indicated at the 
top. We found 1 1 such systems in our sample. 



3 DATA SAMPLE AND OBSERVABLES 

In this section we describe our line identification strategy, absorp- 
tion line measurement techniques and various physically motivated 
subsamples we use for statistical studies. 

For the analysis presented here we concentrate on the inter- 
vening O VI and C iv systems defined as those with apparent ejec- 
tion velocity larger than 50 00 km s~' relative to the QSO emission 
redshift (using Zem given in lRollinde et alj2005l : IScan napieco et al. 
l2006l) . The detailed discussions of t he associa t ed sys tems towards 
QSOs in our sa mple can be found i n Fox et al. liooi). 

Following IScannapieco et al] ( l2006h we define a system by 
grouping together all the components whose separation from their 
nearest neighbor is less than a linking length, Viinjj = 100 km s"'. 
These authors have shown that the absorber's properties and in par- 
ticular the velocity clustering do not change much for velocities 
smaller than this. Note that the same convention was also adopted 
by[Songaila (2005). 

We have searched for O vi following two different approaches. 
(1) First we identify metal line doublets (e.g. Civ and/or Siiv) red- 
shifted beyond the wavelength of the QSO Lya emission and we 
look for O VI doublets around this redshift (i.e. within ~ 100 km s"' 



to the Civ or Siiv redshift). The presence of Ovi is confirmed af- 
ter checking the consistency in the optical depths of the two Ov i 
lines. There are 104 Civ systems (see also D' Odorico et al]|201(]|) 
along the lines of sight we study over the redshift range where O vi 
is detectable with S/N > 10. In addition, there is a system (z^bs 
= 2.7356 towards HE 2347-4342) that is identified by the pres- 
ence of Siiv doublets where the corresponding Civ doublets fall 
in the narrow wavelength range not covered by our UVES spec- 
trum. The presence of Civ in this system has been confirmed by 
lAgafonova et al.l ( l2007h . The total C iv column density in this sys- 
tem is taken from their measurement. Out of these 105 Civ sys- 
tems, 72 show detectable Ovi absorption. Details of the systems 
are given in Table[T] For the other 33 C iv systems, only upper lim- 
its on O VI column density can be obtained. These systems are listed 
in Table |2] Whenever possible we use rms error in the unabsorbed 
continuum at the expected position of the O vi doublets to estimate 
3cr limit on column density using the same number of components 
and fc-values as seen in Civ. In rest of the cases where there is 
strong Lyo- absorption we use the z and fc-values of C iv and gener- 
ated the minimum O vi profile that explains the observed spectrum. 
This allows us to get only a conservative upper limit on ^(O vi). 
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Figure 2. Schematic diagram sliowing the positions of the absorbers in red- 
siiift space along the lines of sight of the QSOs in our sample. The (red) 
circles indicate the redshifts of O vi systems whereas (blue) stars indicate 
the redshift of C iv systems where we have upper limits on ^(O vi). The left 
and right vertical lines are the minimum and maximum redshifts along each 
line of sight. The horizontal dotted lines are just to guide the reader's eye. 



(2) Secondly, we directly search for the O vi doublets in the Lyo- 
forest. In fact, Ovi absorption lines without detectable associated 
C IV and/or Si iv (and H i) can trace the highly ionized gas that orig- 
inates from the WHIM with characteristic temperature ~ 10^-lO^K 
predicted in som e simulations (e.g., Cen & Ostrike3 ll999l. [20o3 : 
iDave et alj[l999ll200lh . For each of these identified coincidences 
we checked the consistency of the shape and optical depth ratios 
of the O VI doublets. We then checked for the presence of associ- 
ated Lyff (and possibly higher Lyman series lines) at the redshift of 
the chosen O vi doublets. While the presence of associated Lyo- ab- 
sorption confirms the O vi identification, it need not be detectable 
in case O vi comes from collisionally ionized gas. Therefore, we do 
not impose the detection of hya absorption as a necessary criterion 
to confirm the Ovi doublets. We found 12 systems from the pres- 
ence of the O VI doublet only (an example is shown in Fig. [TJ. In 
all these cases associated H i is detected. For one of these systems 
(Zabs = 2.7456 towards HE 2347-4342) the wavelength range corre- 
sponding to C IV absorption fall in a spectral gap, so that we could 
not probe the presence of Civ in this system. For the other II sys- 
tems we do not detect any other metal. Apart for the Zahs = 2. 1767 
and 2.3598 systems towards Q 0002-422 and HE 015 1-4326, re- 
spectively, all other systems show absorption from at least one of 
the higher Lyman series lines in addition to Lyo- (i.e. at least Ly/?). 
Thus we are confident that these identifications of O vi are secure. 
However, as we look for the presence of both lines in the doublet 
it is possible that our method has missed some of the O vi only ab- 
sorbers where absorption from one of the transitions (or both) is 
contaminated by Lyor absorption. 

We have found six Lyman-limit systems (LLS) with detectable 
O VI listed in Table[3l There are 3 more LLS (Zabs = 2.7278 & 2.9676 
towards PKS 2126-158 and z.b. = 2.4512 towards HE 0I5I-4326) 
where only upper limits on Ovi column densities could be ob- 



tained. Note that all the LLS are excluded in the analysis presented 
in this paper. This is because A'(H i) measurement in these systems 
are uncertain and they may not trace IGM gas which is of prime 
interest in this study. 

The redshifts of all the intervening Ovi absorbers (both de- 
tections and upper limits) are summarized in Fig.|2l For each line 
of sight (mentioned in the extreme right) the vertical tick mark on 
the left indicates the redshift above which the S/N per pixel of the 
corresponding spectrum is > 10. This defines the minimum red- 
shift (Zmin) for each line of sight as given in column #3 of Table [T] 
& |4] There are 7 systems listed in Table |4] for which the O vi ab- 
sorption falls in the spectral region where S/N < 10. To get robust 
Voigt profile parameters we restrict ourselves to systems detected 
when the spectrum has S/N > 10. The vertical tick mark on the 
right indicates the redshift at which the velocity difference from 
the QSO emission redshift (Zem) is 5000 km s"'. This defines the 
maximum redshift (Zmax) for each line of sight as given in column 
#4 of Table [T]&|4l This cut is applied to remove the Ovi systems 
associated with the QSO or QSO neighborhood. Note that the two 
systems at z^bs = 2.0850 towards HE 1341-1020 and z^bs = 2.1660 
towards PKS 1448-232 falling ~ 5000 km s"' from the respective 
emission redshifts have been included in the sample. 

The maximum redshift path covered by the observations with 
S/N > 10 is Az = 7.62 or AX = 24.85. However, this should be 
treated as upper limits as line blanketing by Lya lines reduce the 
available redshift path length. 



3.1 Absorption line measurement techniques 

3.1.1 Voigt profile fitting 

We use standard Voigt profile fitting and apparent optical depth 
techniques to derive absorption line parameters. The Voigt profile 
fit provides best fitted values of the column density (N), velocity 
dispersion (b) and redshift (z) for each component. The absorption 
lines originating from individual species (Hi, Civ and Ovi) are 
fitted using all the detected transitions with minimum number of 
components required to get the reduced close to I . Whenever 
possible, we have tied the Ovi and Civ components in redshift. 
However, most of the systems are best fitted by components with 
different sets of parameters for O vi and C iv. We use all the avail- 
able Lyman series lines to extract A'CH i). In this case also when- 
ever possible component structure from metal lines were used to 
constrain the redshifts of individual H i components. 

We also have independent Voigt profile decompositions for 
51 Ovi systems along 12 lines of sight (indicated by under- 
lined QSO na mes in c olumn #1 of Table \Q perform ed by 
[Bergeron & Her bert-FortI {2005} using VPFIt]] (WebbI 1 19871 : 
|Rauch et al. 1992). In this case fits to H i, C iv and O vi absorption 
were performed independently without constraining the redshifts 
of any components. As the Ovi absorption lines fall in the Lyo- 
forest, some amount of subjectivity (wavelength range used, place- 
ment of components and number of components etc.,) is involved 
in the Voigt profile decomposition. However, the availability of de- 
compositions derived using two independent procedures allow us 
to investigate the statistical influence of the fitting procedures. 

In total we find 239 individual Voigt profile components for 
Ovi (with 12.75 < log N(Ovi}[cm-^] < 14.49) and 318 compo- 
nents for Civ (with 11.58 < log A'(0 vi)[cm-2] < 14.76). Total 



' See |http://www.ast.cam.ac.uk/~rfc/vpfit.html| 
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Table 4. Details of additional O vi systems detected in the spectral range where continuum S/N < 10 



QSO 








Zsy.s 


log Niii 1) 


log WCOvi) 


log N(C IV) 


Cla,ss 


Case 


<5lypc 


6v{0 VI) 
(kms-') 


(5v(C IV) 
(kms-') 


lAi'lOvi -Civ)| 
(kiTi.s-') 


Low Ions 


Q 0122-380 


2.190 


1.995 


2.137 


1.9746 


15.43+0.16 


14.72+ 0.07 


15.30+0.20 


bd 


B 


-1 


90.4 


281.5 




Yes 


HE 2217-2818 


2414 


1.963 


2.357 


1.9658 


15.82±0.57 


14.95+0.13 


14.42+0.04 


bd 


A/B 





328.3 


171.1 


42.66 


Yes 


Q 0453^23 


2.658 


2.014 


2.597 


2.0041 


15.10+0.09 


14.26+ 0.21 


13.16+0.04 


bd 


B 


+ 1 


166.8 


121.6 




No 


HE 2347^342 


2.871 


2.311 


2.807 


2.1198 


14.26+ 0.04 


14.42+ 0.43 


13.42+0.05 


hb 


A 


+ 1 


34.9 


33.5 




No 










2.2750 


13.97+0.13 


14.59+ 0.77 


13.86+0.20 


bb 


A 


-1 


117.6 


204.0 




No 


PKS 2126-158 


3.280 


2.520 


3.208 


2.4596 


14.25+ 0.05 


14.30+ 0.10 


13.46+0.02 


bd 


B 





48.4 


41.8 


1.55 


No 










2.4855 


15.11+0.16 


14.47+ 0.1 1 


12.61+0.06 


bb 


B 





80.1 


62.2 


16.18 


No 




150 
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3281 328S 3283 3884 3S85 

Wavelength in A 

Figure 3. Bottom : The Ovi /11032 absorption profile of the system at 
Zabs = 2.1808 towards HE 2217-2818. The best fitted Voigt profile is over- 
plotted. Middle : The apparent optical depth profile [t„(v)] estimated from 
the best fitted Voigt profile is shown. The zero velocity corresponds to the 
optical depth weighted redshift vi ~ 2.181378. Top : Integration of the 
apparent optical depth. The line profile velocity width, 6v, has been cal- 
culated as [v(95%) - v(5%)] where v(95%) and v(5%) are the velocities 
corresponding to 95 and 5% percentiles (vertical long dashed lines) of the 
apparent optical depth distribution. 

column densities of Hi, O vi and Civ are listed in Tables [T] and |4] 
in columns #6, #7 and #8, respectively. These are obtained by sum- 
ming the column densities in individual Voigt profile components 
in a given system. 

3.1.2 Apparent optical depth (AOD) 

For each of these systems we have calculated the velocity width 
for O VI an d C iv (wh en detected) absorption using AOD technique. 
Following lLedoux"e t al. ( 2006), the absorption velocity width, Sv, 
has been calculated as [v(95%) - v(5%)], where v(95%) and v(9%) 
are the velocities corresponding to the 95 and 5 percent percentiles 
of the apparent optical depth distri bution. Note thi s definition is 
slightly different from that used by Song ailal ( |2006|) , who defined 
the velocity spread as the velocity range over which the optical 
depth is large r than some frac tion of the peak optical depth. As 
pointed out bv lSongailal J2006t) . the velocity spread measurements 
are very s ensitive to the v elocity range over which the absorption 
is studied. ISongailal l lioo^ used regions within 350 km s ' either 



side of the peak optical depth. This in turn restricts the maximum 
measurable width to 700 km s"' . In order to be consistent with our 
definition of a system, we consider here the velocity range covered 
by the Voigt profile components with no gap larger than 100 km s"' 
(i.e. Viinjj ~ 100 km s"'). 

The (5v(Ovi) is sometimes difficult to measure because of 
blending in the Lycy forest. Therefore, we introduce an index 5iype 
(see column #11 of Tables [T]&|4j which takes the value 0,-1-1 and 
-1. The systems with (5,ype = are the ones for which Ovi ab- 
sorption is well defined based on both the O vi lines. Systems with 
Stypc = + 1 and - 1 are those for which the measured width should 
be considered as an upper and lower limit, respectively. There are 
54 systems with Siypc = 0. 

We also calculate the first moment of the apparent optical 
depth distribution which gives us the optical depth weighted red- 
shift for both Civ and Ovi (i.e. zqi^ and zqvi)- This allows us 
to calculate the velocity shift (called |Av(0 vi - C iv)|) between the 
O VI and C iv absorption originating from the same system. 

Fig.[3]illustrates the way we measure dv. Bottom panel shows 
the observed Ovi profile and the best fitted Voigt profile. In the 
middle panel we plot the AOD profile as a function of velocity. 
The top panel shows the integrated optical depth starting from a 
minimum velocity. The left and right vertical dashed lines mark the 
velocities at which the integrated optical depth is 5% and 95% of 
the total optical depth respectively. The velocity difference between 
these two lines gives 6v. For clarity, the zero velocity is fixed at 
zq so that integrated optical depth is 50% at v = km s"' . 

If both O VI and C iv absorption originate from the same gas 
with constant density and homogeneous ionization conditions then 
the velocity offset measured between zqvi ^Civ should be 
zero (i.e. |Av(0 vi - Civ)| = 0). Any mismatch of the optical depth 
weighted redshifts between species could be interpreted as (a) rel- 
ative line of sight velocity between the two species (as expected in 
the case of O vi originating from different interfaces like evaporat- 
ing region, cooling front or shocked gas) or (b) ionization inhomo- 
geneity along the line of sight. The measured Sv(0 vi), 6v{C iv) and 
|Av(0 VI - C iv)| in our sample are summarized in column #12, #13 
and #14 of Tables [T]&|4] respectively. 

3.2 Sub-samples of O vi systems 

As O VI originates from a whole range of physical conditions and 
because of blending, the characteristics of the O vi absorption are 
difficult to extract. Hence, we build different subsamples depending 
on the issue we want to address. 

The H I column density and fe-parameter are very important 
for un derstanding the physical state of O vi gas (see e.g., Tripp et al] 
|2008|) . In particular, to investigate whether the gas is photoionized, 
we need to know A'(H i) in individual O vi components. So we de- 
fine "Case- A" systems as the ones where at least one of the Lyman 
series line is unsaturated. In principle these systems also enable us 
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Figure 4. Examples of O vi systems where both Unes in the doublet are unblended (left; dd), one of the Hne in the doublet is blended (middle; bd) and both 
lines in the doublet are partially blended (right; hh). The smooth curves are the best fitted Voigt profiles. The horizontal tick marks indicate the centroids of 
the individual components. The absorption redshift that defines the zero velocity and the name of the background QSO are indicated at the top of each panel. 



to derive Z?-values for most of the Hi components. There are 35 
systems classified as "Case-A" systems in our sample. The systems 
shown in the left and middle panels of Fig.|4]are examples of "Case- 
A" systems. There are 10 systems where at least one Hi absorption 
is unsaturated for some of the components. We call these systems 
at "Case-A/B". The right most panel in Fig.|4]is an example of such 
a case. In the remaining systems most of the H i absorption is sat- 
urated and Voigt profile measurements are uncertain. We call these 
systems "Case-B". We notice that several of the "Case-B" systems 
show absorption lines from low ionization species such as C ii, C iii, 
Si m/Si IV etc. The information regarding the presence of low ions is 
also given in column #15 of Tables[T]|2]and|4] In total there are 46 
systems where low ions are detected. Note that the C iii information 
is not available in our sample for z ^ 2.2. 

We also classify the systems based on the presence of un- 
blended Ovi doublets. Out of 84 systems only 29 systems have 
both O VI lines unblended for which we are able to perform robust 
Voigt profile fits. We call this sample as "robust sample" as far as 
Voigt profile parameters are concerned (these are marked as "rfd" 
in column #9 of Tables [T] and |4l(. An example of such a system is 
shown in the left panel of Fig.|4l In 30 other systems one of the O vi 
lines is blended (as in the middle panel of Fig.O. Here the O vi fit 
relies mainly on the unblended line. We denote these systems as 
"W in column #9 of Tables [T] and O In the remaining 25 systems 
both lines are partially blended (as in the right panel of Fig.O. In 
this case (referred to as ''bV in column #9 of Tables [T] and |4l( we 
use part of the profiles of both O vi lines to fit the Voigt profiles. 
In most of these cases the additional absorption is treated as due to 
intervening Lya absorption. Therefore, the Voigt profile parameters 
(number of components and fc- values) are uncertain. 



4 ANALYSIS BASED ON VOIGT PROFILE FIT 

In this section we discuss the distributions of various parameters 
derived from our Voigt profile fits. 



4.1 ii-parameter distribution 

The Z7-parameter derived from Voigt profile fit provides only an up- 
per limit to the kinetic temperature of the gas. Thus, we do not 
attempt to constrain the physical state of the O vi gas using indi- 
vidual ii- values, instead we draw some broad conclusions using the 
/7-parameter distributions. 

First, we compare our O vi fe-parameter distribution for t he 12 
lines of sight with that of the Bergeron & Herbert-Fort (2005). The 
Kolmogorov-Smirnov (KS) test suggests that the two distributions 
are drawn from the same parent population with a 32% probability 
for the observed deviation in the cumulative distributions to occur 
by chance. 

In the bottom panel of Fig. |5] we plot the fc-parameter distri- 
butions measured in individual Voigt profile components for O vi 
(solid histogram) and C iv (dot-dashed histogram) in our full sam- 
ple. The (red) dashed histogram gives the ^1(0 vi) distribution in the 
robust components (i.e. from "dd" subsample). The Ovi compo- 
nents arising from systems with S/N < 10 (listed in Table|4l( are not 
included here. Median values of b(0 vi) = 13.8 km s"' and b{C iv) 
= 10.1 km S-' correspond to T ~ 1.8 x 10' K and 7.4 x 10* K for 
O VI and C iv respectively in the case of pure thermal broadening. 
The temperature corresp onding to median b{0 vi) agrees well with 
r ~ 2. 1 X 10' K found bv lSimcoe etakl ( l2002h . The median fc-value 
of O VI components becomes 12.5 km s"' if we restrict ourselves to 
the robust sample. In the case of pure photoionization heating the 
expected temperature (i.e. T ~ 2 x 10* K) implies fe(Ovi) ~ 5 
km s"'. The lowest b(0 vi) measured (4.9±3.2 km s"') in our sam- 
ple is consistent with this. 

In the case of collisional ionization, the fraction of oxygen 
in Ovi peaks around T ~ (2 - 3) x 10' K, which corresponds to 
fc(Ovi) ~ 14.4 km s"'. This is very close to the median fc(Ovi)- 
value of our full sample. This is shown by the left vertical dashed 
line in Fig. [5] There are 52% (and 62% for the dd subsample) com- 
ponents having ^-parameters less than 14.4 km s"' suggesting pho- 
toionization (and/or non-equilibrium collisional ionization at high 
metallicities) is the dominant process in these systems. On the other 
han d, the O vi fraction will b e less than 0.01 when T > 6 x 10' K 
(see lGnat& Sternberg! l2007h . or fe(Ovi) > 25.5 km s ' (second 
vertical line in Fig.O. Therefore under the purely thermally broad- 
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b-parameler [km s"'] b(0 VI) [km s"'] 



Figure 5. Bottom : Comparison of the fo-parameter distributions of C iv (dot- 
dashed histogram) and O vi (solid histogram). The dashed histogram is for 
the O VI components originating from dd systems. The vertical dashed lines 
correspond to b(0 vi) = 14.4 and 25.5 km s"' respectively. Top : The cumu- 
lative /^-parameter distributions of C iv, O vi (both total & dd) components 
are shown following the same line style as in the bottom panel. 



ened case one does not expect the O vi fc-parameter to be higher 
than this unless the metallicity and/or N(H) is very high. These 
systems will also have broad and shallow associated Lyor. We find 
14% (and 8% for the dd subsample) of the components having b- 
values higher than this. Most of these high b components in our 
sample are part of blends where the O vi profile is decomposed into 
multiple Voigt profile components. Isolated O vi components with 
b{0 vi) > 25 km s"' are very rare. Note that very few such isolated 
broad O vi components together with broad albeit shallow Lyg ab - 
sorption are detected at low-z (see e.g.. , Savage et al..,201oll2oTlh . 
These are interpreted as collisionally ionized gas with T ~ 10* K. 
The non-detection of such systems in our sample may be attributed 
to the bias introduced by line blanketing and blending due to the 
Lyo- forest absorption. 

The temperature for which the ionization fraction of Civ 
peaks under collisional ionization equilibrium is ~ 1.1 x 10^ K, or 
biC iv) = 12.40 km s"'. Sixty seven percent of the C iv components 
have fc-parameter less than this suggesting that a considerable frac- 
tion of C IV originates either from photoionized gas or from non- 
equilibrium cooling gas with high metallicity (i.e. Z ~ 1.0 Zq). It 
is apparent from Fig. |5] that the fo-parameter distribution is wider 
in case of O vi compared to C iv. This is evident in the cumulative 
distributions plotted in the top panel. The KS-test indicates that the 
two distributions are drawn from significantly different populations 
with the probability of this difference occurring by chance being 
less than 0.1%. Ideally one would expect fc(Civ) > fo(Ovi) if Civ 
and O VI trace the same gaseous phase. Our finding is consistent 
with O VI and C iv originating from different p hases of the gas as- 
sociated with the absorption sy stem (see also, Isimcoe et al I2OO2I : 
[Bergeron & Herbert-Forll20051) . 

Next we compa re the fc(Ov i ) distr ibution in our sample with 
the low-z sample of iTripp et al. I ( l2008h obtained using the Space 



Figure 6. Comparison of the /^-parameter distiibutions of O vi at high and 
low re dshifts. The low-z sample is the STIS data sample of iTripp et al J 
Only those Ovi components, originating from systems with 
log ^(0 vi)(cm"^) > 13.6 (corresponding to a rest frame EW of ~ 50 mA) 
are considered here. The vertical dashed hnes are as in Fig. [5] The overall 
/:)-distribution of our sample is significantly dift'erent from that of the low-z 
sample as can be seen from the cumulative distributions plotted in the top 
panel with the same line style as in the bottom panel. 



Telescope Imaging Spectrograph (STIS) on board HST. The spec- 
tral resolution of STIS data is ~ 7 km s"', very much similar to 
the one used here but the STIS spectra are slightly under-sampled. 
Even though the S/N of STIS spectra are typically a factor of two 
lower than those of our UVES data, the contamination by interven- 
ing H I absorption is less severe in the low-z sample. In the bottom 
panel of Fig. [6] we show the two O vi ^-distributions. We restrict 
the comparison to systems with log A'(Ovi) (cm^^) > 13.6 which 
is the completeness limit of the STIS sample. The median Z?- values 
for the low and high-z Ovi components are 24.0 and 14.6 km s"' 
(13.2 km s"' for the dd subsample), respectively. It is interesting to 
note that almost 87% of the low redshift fc-parameters are consis- 
tent with the temperature expected from collisional ionization, i.e., 
i(Ovi) > 14.4 km s"' (first dashed vertical line), only 53% (44% 
for the dd subsample) of the high redshift components satisfy this. 

Almost 43% of the low redshift components are consistent 
with b > 25.5 km s"' (second dashed vertical line). In our high 
redshift sample, only 16% (9% for the dd subsample) components 
show fc- value greater than 25.5 km s"' . The overall ^-distribution of 
our sample is significantly dilferent from that of the low-z sample 
as can be seen from the cumulative distributions plotted in the top 
panel of the Fig. |6l A two sided KS-test gives a maximum depar- 
ture between the two distributions, D = 0.41 with a probability, P 
= 3.9 X 10"* that the two samples are drawn from the same parent 
population. The difference is even more when we compare the low- 
z sample to our dd subsample as can be seen from the cumulative 
distributions. 

In summary, we find that at high redshift, the fo-parameter 
distributions of Civ and Ovi are significantly different with Ovi 
absorption being wider than C iv, suggesting that the two species 



High redshift O vi absorbers 9 



Table 5. Results of decomposition of thermal and non-thermal contributions to the line broadening 



QSO 


Zsys 


Vrcl 




h(0 VI) 


b(C IV) 


log 


; N (era 






logr 


log Tniax'^ 






(km s-') 


(km s-') 


(km s-') 


(km s-') 


Hi 


Ovi 


Civ 


(km s-') 


(T in K) 


(r,„ax in K) 


PKS 1448-232 


2.1099 


-4.7 


14.2 


8.1 


7.3 


13.05 


14.28 


13.11 


7.5 (6.3) 


3.94 (3.99) 


<4.09 


PKS 0237-23 


1.9878' 


-hO.4 


20.6 


9.5 




13.64 


13.41 




8.2 


4.33 


<4.4I 


PKS 0237-23 


2.0108' 


-3.8 


31.3 


9.4 




14.31 


13.15 




5.4 


4.76 


<4.77 


Q 0329-385 


2.0764 


+4.9 


19.7 


7.9 


6.2 


13.64 


13.26 


13.21 


6.4 (2.6) 


4.32 (4.36) 


<4.37 


Q 0329-385 


2.2489' 


-50.3 


18.6 


8.2 




12.76 


13.68 




7.0 


4.26 


<4.32 


Q 0329-385 


2.2489' 


-19.4 


10.1 


10.3 




11.95 


13.47 








<3.79 


Q 0329-385 


2.2489' 


- 3.1 


33.1 


21.9 




13.11 


13.76 




20.9 


4.60 


<4.82 


Q 0329-385 


2.2489' 


+43.4 


24.7 


12.1 




13.10 


13.64 




10.7 


4.48 


<4.57 


Q 0329-385 


2.3139 


+4.0 


38.6 


11.5 


4.3 


14.02 


13.16 


11.79 


6.5 (...) 


4.94 (...) 


<4.96 


Q 0329-385 


2.3139 


+34.8 


16.9 


8.8 


5.8 


13.23 


12.75 


12.21 


8.0(3.3) 


4.13 (4.22) 


<4.24 


Q 0329-385 


2.3639 


-46.3 


27.8 


15.3 


9.2 


14.40 


13.37 


11.86 


14.1 (4.7) 


4.54 (4.66) 


<4.67 


Q 0329-385 


2.3639 


-4.1 


19.0 


8.7 


9.1 


13.91 


13.52 


12.25 


7.5 (7.6) 


4.27 (4.26) 


<4.34 


HE 1347-2457 


2.3327' 


+3.2 


23.5 


18.5 




13.78 


13.26 




18.1 


4.13 


<4.52 


HE 1347-2457 


2.3422' 


+ 10.8 


18.3 


9.9 




14.40 


12.95 




9.1 


4.18 


<4.31 


Q 0453^23 


2.5371' 


-26.3 


38.7 


10.3 




13.48 


13.05 




3.6 


4.95 


<4.96 


Q 0453-423 


2.5371' 


-1.6 


28.5 


11.7 




14.62 


13.17 




9.6 


4.64 


<4.69 


PKS 0329-255 


2.5687 


-3.9 


30.0 


10.1 


6.1 


14.62 


13.12 


12.15 


7.0(...) 


4.71 (...) 


<4.74 


PKS 0329-255 


2.5687 


+ 14.2 


29.9 


13.8 


20.8 


14.25 


13.35 


12.27 


12.0(19.8) 


4.66 (4.48) 


<4.73 


HE 0151-4326 


2.5053 


+42.6 


28.0 


21.7 


10.5 


14.71 


13.52 


12.10 


21.2 (7.0) 


4.31 (4.65) 


<4.68 


HE 2217-2818 


2.0748 


+64.4 


29.0 


12.6 




14.12 


14.32 




10.6 


4.64 


<4.71 



Table Notes - ' Components from "O vi only" systems. ^Calculated from b(H i) assuming pure thermal broadening. 
Values in the parenthesis are calculated using h(}i i) - h(C iv) pairs. 



trace different phases of the absorbing gas. The fc(Ovi) distribu- 
tio n at high-z i s very different f rom t hat at low-z as measured 
by iTripD et alj ( |2008|) . Recently |Fo3 ( |2011|) has drawn a sim- 
ilar conclusion using th e high -z Voigt profile fitting results of 
iBergeron & Herbert-Fori ( l2005l) . 



4.2 Thermal and non-thermal contributions to ^-parameters 

In the previous section we find that on an average fc-parameters 
of O VI measured at low-z are higher than that measured at high-z. 
To explore this further we use a subsample of O vi absorbers where 
O VI and H i absorption are well aligned. In general the line width of 
any species can be decomposed into thermal (fc,h) and non-thermal 
(but) parts, i.e. = + b^^. For species located in the same 
physical region, the non-thermal part is supposed to be identical 
whereas the thermal part scales inversely with the mass of the ion. 

Line saturation in the case of Lyo- and blending with Hi lines 
in the case of O vi make the robust estimation of ^-parameters for 
Hi and Ovi difficult in our sample. In addition we need to ensure 
that there is a good alignment between Ovi and Hi absorption. 
Hence we select systems using the following two criteria : (1) The 
component structure of Hi is well defined and one of the avail- 
able Lyman series lines is unsaturated (i.e. "Class-A" absorbers as 
defined in Table [T). (2) The Ovi profiles are well defined and the 
velocity offset between H i and O vi absorption centroids are con- 
sistent with zero within 3cr uncertainty. We also avoid systems with 
low ion absorption lines as Hi seems to be predominantly associ- 
ated with the low-ionization phase when these species are present 
(see section|6j. Thus by using these selection criteria we minimize 
the probability that these absorbers have a multiphase structure. 

We find only 13 systems with 19 Voigt profile components 
(identified by vertical dashed lines in Fig. |7]l satisfying the condi- 
tions listed above. This is only ~ 15% (i.e. 13 out of 84 systems) of 
our full sample. There are 6 systems showing only O vi and H i ab- 
sorption (i.e. "O VI only" systems). In the remaining 7 systems both 
O VI and C iv absorption are seen. In 6 of them O vi and C iv compo- 



nents are remarkably aligned. In these cases we estimate the tem- 
perature and b„t using both fc(H i)-b(0 vi) and i(H i)-b(C iv) pairs. 
Only in the case of the Zab.s= 2.0748 system towards HE 2217-2818, 
the corresponding Civ component is ~ 4.2 km s"' away from the 
O VI component and hence it is not used in our analysis. 

The results of the decomposition in thermal and non-thermal 
broadening are summarized in Table|5] Columns #1, #2 and #3 list, 
respectively, the QSO name, system redshift (Zsys) and the veloc- 
ity of the component (vrei) with respect to the systemic redshift. 
Columns #4, #5 and #6 give ii-values of H i, O vi and C iv compo- 
nents, respectively. The corresponding column densities are given 
in columns #7, #8 and #9 respectively. The non-thermal contribu- 
tion to the broadening and the estimated temperature are listed in 
columns #10 and #11, respectively. The values in parenthesis are 
calculated using fe(H i)-ii(C iv) pairs. Column #12 lists the upper 
limits on the temperature as calculated from fe(H i) assuming pure 
thermal broadening. It can be seen from the table that the temper- 
atures estimated from the b(lii)-b{Ci\) pairs are consistent with 
those derived from the fc(Hi)-fo(Ovi) pairs, whereas fo(Civ) gives 
slightly lower fc„t compared to that obtained using ^(O vi). 

In Fig. [8] we show the distribution of (top) and temperature 
(bottom) as calculated from the fc(H i)-b{0 vi) p airs. The solid his- 
togram shows the results bv lTripp et ai] l l20081) for the low-z well 
aligned O vi absorbers for comparison. The values of non-thermal 
velocity in our sample are found to be in the range 3.6 < fcm ^ 21.2 
km s"' with a median valu e of 8.2 km s"' . The median value of b„t 
for the iTripD et al.l ( |2008|) sample is ~ 20.0 km s ' . From the top 
panel of Fig. [8] it is apparent that the b„i distributions at high and 
low-z are significantly different. This is confirmed by a two sided 
KS test with a probability that the two distributions differ much 
greater than 99.9%. 

The median value of the temperature distribution in both high 
and low redshift samples is found to be ~ 3 x 10"* K. The KS 
test shows that the temperature distribution of low-z sample is not 
significantly different from high-z sample (~ 38% probability that 
the difference is occuring by chance). It is interesting to note that 
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Figure 7. Velocity plots of the 13 well aligned systems that are used to measure temperature and turbulent motion of the absorbing gas. The best fitted 
Voigt profiles are over-plotted on the observed data. The vertical tick marks show the positions of individual components. The vertical dashed lines mark 
the components used for this study. The occasional dashed curves over-plotted on the data show the contamination from other intervening absorption. The 
absorption redshift that defines the zero velocity and the name of the background QSO are indicated at the top of each panel. 
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Table 7. Summary of power law index measurements for C iv CDDF 

Redshift > log N(C iv) pQ jy References 

Systems Components 




a 

s r 

o - ' 

^ 6 - I I 

3 - ' 

1 ■ ; ■ 

^ L- 1 - -I - J I I I 1 I 1 I ! I 1. - i - X - iZl- 

3.6 3.9 4.2 4.5 4.8 5.1 5.4 5.7 
log T (in K) 

Figure 8. Distributions of temperature (lower panel) and turbulent velocity 
(upper panel) in the well aligned sample are shown in dashed histogram. 
The solid histogra m in both panels sh ows the corresponding distribution at 
low-z as derived bv lTripp et aLl i2008l) . 



1.9-3.1 12.6 1.6±0.1 1.9±0.1 1 

1.9-3.1 13.0 1.7±0.1 2.0±0.1 1 

13.4 1.64±0.10 .... 2 

2.9-3.5 13.0 1.8±0.1 .... 3 

13.0 1.71±0.07 .... 4 

12.3 1.44±0.05 .... 5 

1.6 1.8 6 

< 1.0 13.2 1.50+!J!' ... 7 



References:- (1) This work; (2) Petitjean & Bergeron 1 1994); (3) 
Songaila (2001); (4) D'Odorico et al. (201 01): (5)|Ellison et'alj i2000l) : (6) 
iBoksenberg et al.. 1.2003.) ; (7) .Cooksev et alH201ol) 

high correlation coefficient (> 0.8) and at > 3(T significance level 
in both samples. A weak correlation (at < 2a level) is seen be- 
tween b^f a nd A^(0 vi), both at hi gh and low redshift. Recent sim- 
ulations of ICen & Chisaril ( 1201 ll) suggest such a trend of higher 
non-thermal contribut i on at higher A'(0 vi). On the other hand 
lOppenheimer &Davd ( l2009l) have introduced density dependent 
turbulence in their simulations in order to reproduce the equivalent 
width distribution and the b-N correlation of low-z O vi absorbers. 
The lack of strong correlation between and A'(0 vi) or A'(H i) in 
our sample suggesting that this density dependence of turbulence 
may be weak at high-z. 



Table 6. Results of Spearman rank correlation analysis 



Observable- 1 


Observable-2 


This Work 


TrioD et al. (2008) 


P.S 


pjo- 


Ps 




i'nc(Ovi) 


^(O VI) 


+0.41 


+ 1.7 


+0.31 


+ 1.6 


font (0 VI) 


A'(Hi) 


+0.21 


+0.9 


+0.07 


+0.4 


fcnl(Ovi) 


b{0 VI) 


+0.78 


+3.3 


+0.99 


+5.2 


fcnl(Ovi) 


fc(Hi) 


-0.05 


-0.2 


+0.56 


+2.9 


log 7(0 VI) 


fc(Ovi) 


+0.30 


+ 1.3 


+0.21 


+ 1.1 


log 7(0 VI) 


fo(Hi) 


+0.92 


+3.9 


+0.83 


+4.3 



while none of the components has temperature log T > 5.3 (which 
would favor collisional ionization for O vi), 42% of the components 
(i.e. 8 out of 19 components) show 4.6 < log T < 5.0, which is 
warmer than the temperatures expected in photoionization equi- 
librium. These higher temperatures can be obtained in a rapidly 
cooling over-ionized gas that was shock heated through mechan- 
ical processes such as galactic winds. We thus compare the ob- 
served A'(0 vi)/A'(Hi) in these well aligned comp onents with the 
non-eq uilibrium collisional ionization models of Gn at & Sternberg 
( I2OO7I) assuming the temperature derived from b(ll i) (given in col- 
umn #12 of Table O and find that the observed ratios can not be 
reproduced by these models (even when we use maximum gas tem- 
perature) for gas phase metallicity less than solar. Thus it seems that 
the ionization state is probably maintained by the UV background 
radi ation which is expecte d to dominate when T < 10' K (see Fig. 
9 of iMuzahid etalj201ll) . 

In Table|6]we summarize the Spearman rank correlation anal- 
ysis to search for possible correlations between log T and 
other observables of Ovi absorbers discussed here as well as in 
iTripD et al] ( I2OO8I) . In all the correlation analysis presented in this 
paper, the correlation coefficient and its significance are denoted 
by Pj and pjcr, respectively. As expected, b„t and log T are 
strongly correlated to b(0 vi) and Z?(H i), respectively, with a very 



4.3 Column density distributions 

In this section, we study the column density distribution functions 
(CDDF) of Civ and Ovi systems in our sample. We use the usual 
parametrization of the column density distribution function, i.e., 

f(N)dN = BN-l^dN, (1) 

where, fiN) is the number of systems/components per unit column 
density interval per unit redshift path length defined as, 

dX = (i+ zf[Q.K + £lm(l + zfr"^dz. (2) 

We use the maximum likelihood method to estimate the power 
law index, p. The column density distribution functions of O vi 
systems (top) and components (bottom) are shown in the left 
panel of Fig. |9] The (blue) squares in the bottom panel corre- 
spond to the Ovi component CDDF derived from the fits by 
[Bergeron & Herbert-FortI ( |2005|) for only 12 sight lines. It is in 
good agreement with our results based on the full sample. As men- 
tioned above, we cover a redshift path (Az) of 7.62 (and /S.X = 
24.85) with the full sample where the S/N is > 10 per pixel. This 
figure shows that our survey is not severely affected by incomplete- 
ness for log ^(O vi) > 13.7 as indicated by the vertical dotted line. 

A maximum likelihood fit to our data gives a power law index 
of ySQ VI = 1 -9 ± 0. 1 for O VI systems and pQ = 2.4 ±0.2 for O vi 
components for log A'(0 vi) > 13.7 . We would like to mention that 
in the subsample dd where both lines in the doublet are unblended, 
we find /Jqvi = 2.2 + 0.3 for systems and jSovi ~ 2-7 ± 0.4 for 
components . These values are somewhat steeper than, albeit con- 
sistent with, the full sample within the measurement uncertainties. 
For a sample of low redshift (z < 0. 15) O vi absorbers with Lya res t 
frame equivalent width (VKLy„) > 80 mA. lDanforth & Shulll j2005l) 
found Pq yj = 2.2 ± 0. 1 for the componen ts. In a latter paper study- 
ing an extended sample of O vi absorbers, iDanforth & ShullH2008l) 
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Figure 9. Left: Column density distribution functions of vi systems (upper panel) and Ovi components (lower panel). The y-axis is the number of sys- 
tems/components per column density interval per unit redshift path length. The bi n size along the column density axis is 10''^cm"^ and the Icr error bai's in 
the y-axis are calculated using Poisson statistics. The (blue) squares are taken from lBergeron & Herbert-Forll )200S) . The dashed straight lines are power laws 
of the form f(N) = BN^^. The power law indices, /3, are obtained using the maximum likelihood method where no binning is involved. Vertical dotted lines 
show the column density above which our sample is complete. Right: Same as left but for C iv. 



found jSqvi to be 1.98±0. 11. In addition, we have performed a 
maxim um likelihood analysis on the low-z sample of Tripp et al. 
( I2OO8I) and found ySQ^j = 2.3 ± 0.2 for components and 2.0±0.2 
for systems for log A'(Ovi) > 13.7. Therefore the jSqvi measure- 
ments at low-z are consistent with our estimations. 

The right panels of Fig.|9]show the C iv column density distri- 
bution function for components (lower panel) and systems (upper 
panel). Incompleteness limit is clearly log A'(Civ)(cm"^) > 12.6 
as shown in the figure by a vertical dotted line. Using maximum 
likelihood analysis, we find, Pqi^, to be 1.6±0.1 for systems and 
1.9±0. 1 for components for log A^(C iv) > 12.6. Our results are 
consistent with the earlier findings jPetitiean & BergeronI 1 19941 : 
ISongailall200ll : lD'Odorico et alj201(]h within Icr uncertainty when 
we use similar column density cutoff (i.e. log A'(Civ) > 13.0; see 
Table|7J. 



4.4 O VI and C iv cosmological densities 

In this section we compute the contribution of Ovi and Civ ab- 
sorbers to the baryon density using our column density estimates. 
The cosmic density of the O vi absorbers can be expressed as : 



Ovi 



WomoUSiMCOvi) 



c p,., 



AX 



(3) 



with an associated fract ional variance (as proposed by 
IStorrie-Lombardi et al.lll996h : 



^OyJ [2.A'i(0vi)p 



(4) 



where Hg is the Hubble constant, mo is the atomic mass of oxygen, 
AX and pc, are the total redshift path and the current critical density, 
respectively. We find IIq^j = (1.0 ± 0.2) x 10"^ for log A'(0 vi) > 
13.7. If we include the LLS in our calculation then f^Qyi = (1-3 - 
0.3) X 10"' for the same column density cutoff. Usi ng the column 
densities of 12 intervening systems listed in Table 2 of lSimcoe et alj 
f2002) we calculate Q.q yj for their sample. For their uncorrected 
redshift path (i.e. AX = 6.9), I^Qvi out to be (1.2±0.5)x 10"'' 
for logA'(Ovi) > 13.7 which is in good agreement with what we 
find here. For the 12 lines of sight our estimated value of fl p is 
in excellent agreement with 'Bergeron & Herbert-For3 ( l2005h . 

It is also possible to calculate the fractional contribution to the 
cosmological density of the baryons associated with the O vi phase, 
provided the ionization fraction of Ovi (/qvi) ^id the metallicity 
(Z) of the gas are known, viz.. 



Ovi ^(Hoj^m] 



tPcr A,/OviZ(0/H)oA AX 



2iA'i(Ovi) 



(5) 



Where fj = 1.3 is the mean atomic weight and nin is the mass 
of the hydrogen atom. For log A'(Ovi) > 13.7 the value of fiS,^' 



Ovi 



0.001 1 [(/V0.71)(Z/0. lZo)(/o vi/0.2)]-' . For /q , 



is : Q. 

0.2 and Z = O.IZq, we get flOjVi^b = 0.028 i.e. Ovi absorbers 
contribute 2.8% to the total baryon density at redshift z ~ 2.3. For 
a more typical IGM metallicity of Z = O. OlZg this contribution can 
go up to ~ 30%. At low redshift (z < 0.4)' Danforth & Shulll ( |2008|) 
derived 0.2^/^^ = 0-073 ± 0.008 down to log N{0 vi) = 13.4 and 
0.086 ± 0.008 down to log A'(Ovi) = 13.0 (for /q^j = 0.2 and 
Z = O.lZo). 

Using Eqn. {Sjl for log ^(C iv) > 12.6, we find CIq = (2.4 ± 
0.6) X 10"** without including Lyman limit systems. Inclusion of 
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Table 8. Summary of CIq jy measurements 



Table 9. Results of correlation analysis between ^(O vi) and b{0 vi) 



Redshift 


nc,vX(10-») 


References 


< 1.0 


6±1 


1 


1.9-3.1 


5.3±2.2 


2 


2.0-2.5 


6.83 


3 


2.0-2.5 


5.0±1.0 


4 


~ 2.2 


7.5±2.2 


5 


~ 2.5 


3.8±0.7 


6 


2.5-3.0 


4.79 


3 


2.5-3.0 


3.2±0.7 


4 


5.3-6.0 


< 1.3 


7 


~ 5.8 


0.44±0.26 


8 



(2) This work; ( 3 ) ISongailj 



(2001.); (4) .D'Odorico et al, 12010.); (5) 


Scannapieco et al.l J2006l); (6) 


iBoksenbers et al.l J2003h; (7) iBecker et al 


J2009l); (S) iRvan-Weber et alj 
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Figure 10. ^(O vi) against b(0 vi) measured in individual components. The 
short dashed curve and the dotted line show the So" detection threshold and 
the N(0 vi) limit for the sample completeness respectively. 



LLS increases the value up to (5.3±2. 1)xl0 ^ . This value compares 
well with the previous estimations ( Songaila ' 200ll; Isimcoe et aP 
l2004IScannaDieco et aI.l2006l : lD'Odorico et alT2010t) for a similar 
redshift range as can be seen from Table (8] For completeness, the 
recent SIq jy measurements for z > 5 are also summarized in this 
table. 



4.5 b - N correlation 

A correlation between column density and ^-parameter of Ovi 
compo nents has been predicted by theoretical studies of hot ionized 
gas (see lEdgar & Chevalieill98^ : lHeckman et alj2002l) as a natural 
consequence of radiatively cooling hot gas passing through a coro- 
nal regime. Heckman et al. ( 2002) have shown that such a correla- 
tion exists in a wide variety of astrophysical environments (such as 
Galactic disk, halo. High Velocity Cloud (HVC), Large Magellanic 
Cloud (LMC), Small Magellanic Cloud (SMC), starburst galaxies, 
IGM etc.). These authors have also shown that the relationship be- 



Sample 


> log A'(Ovi) 


Components 




Ps/tr 


Total 


13.7 


82 


0.05 


0.4 


dd 


13.7 


19 


-0.36 


-1.5 


bd 


13.7 


22 


0.17 


0.8 


bb 


13.7 


41 


0.24 


1.5 



tween the log ^(O vi) and log b is linear for the broad lines (i.e. 
fc(Ov i) > 40 km s~' ) but ro l ls over and steepe ns for the narrower 
lines. iDanforth et al.l JioO^ . lTripp et"al] ( |2008[) on the other hand, 
report no convincing evidence of such correlation for the low red- 
shift intergalact i c O vi absorbers. iLehner et al.l ( [2006h revisited the 
iHeckman erZI j2002l) model and found that their sample of O vi 
absorbers is consistent with the model but the observed A^(Ne viii) 
is much less than the model prediction. 

In Fig. [To] we plot the Ovi ii-values against column densities 
in individual components for our full sample. The dashed curve 
shows the 5cr detection threshold of our data assuming S/N ~ 10 in 
the forest. The limiting equivalent width f or a given fe-value ha s 
been calculated using the prescription by iHellsten et al.l ( Il998h . 
This limiting equivalent width is then converted to a column den- 
sity assuming the optically thin case. The vertical dotted line at 
log A'(Ovi) = 13.7 shows our sample completeness (see section 
|4.3l l. For log A'(0 vi) > 13.7 the b-N space is uniformly populated 
by the data points indicating the lack of any significant correla- 
tion. The results of Spearman rank correlation analysis performed 
between N(0 vi) and b{0 vi) in our full sample and various subsam- 
ples are given in Table |9] No statistically significant correlation is 
found in any of these cases when appropriate column density limit 
(i.e. log A'(0 vi) > 13.7) is considered. 



5 ANALYSIS BASED ON APPARENT OPTICAL DEPTH 

In this section, we study the line kinematics of O vi and C iv ab- 
sorption using the line spread (6v) and the velocity shift (Av) be- 
tween the C IV and O vi optical depth weighted redshifts defined in 
section [3. 1.21 For most of the discussions presented below, unless 



otherwise stated, we restrict ourselves to the S, 
defined in Table[T] 



type 



subsample 



5.1 Line spread (Sv) distribution 

In the left panel of Fig.[TT] we show the <5v distribution of O vi sys- 
tems (solid histogram) and C iv systems (dashed histogram, includ- 
ing all C IV systems). The median value of Sv(0 vi) is 66 km s"' and 
the maximum observed value is 340 km s"'. The median value of 
^i'(Civ) is 58 km s"'. The Sv distributions of Ovi and Civ, ap- 
pear fairly similar except in the first velocity bin. If the latter is 
considered, the KS-test gives a probability of 0.1% that the two 
distributions are drawn from the same population (see top panel in 
Fig. II lb. However if it is excluded, then the KS-test shows only a 
~ 52% probability that the two distributions differ. The lack of O vi 
systems with dv(0 vi) < 40 km s"' could be related to the dilficulty 
in detecting narrow systems in the Lya forest. 

In the right panel of Fig. [TT] we plot the C iv line spread 
against that of Ovi. The (red) filled circles are for systems with 
^type = 0. A strong correlation between 6v(C iv) and ^v(0 vi) is ap- 
parent from the figure. The Spearman rank correlation coefficient 
for the robust measurements is ~ 0.81 with ~ 5.3cr significance. 
The slope (0.8 1 ± 0. 1 1) and the intercept (-8.60 ±9.17) of the best 
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Figure 11. Left : Sv distributions of Ovi (solid liistogram) and Civ (daslied liistogram) systems. Corresponding cumulative distributions are plotted in the 
upper panel. Right : (5i'(Civ) versus (5i'(Ovi). The (red) solid circles are for the systems with i5|ype = whereas the an'ows indicate systems with (5typc = ±1. 
The straight line shows the least square fit to the data for systems with (5|ype = 0. 



fitted straight line to the (red) solid circles indicates that the spread 
of C IV absorption is systematically smaller than that of O vi for a 
given (5v(0 vi). 

There are three O vi (top panel of Fig. [T2l ( and three C iv sys- 
tems (bottom panel of Fig.[T2j with 6v > 200 km s"'. These large 
velocity spreads could be related to either (a) large scale wind s 
as seen in Lyman break galaxies l lAdelb erger et al. 2003, 200?), 
(b) redshift clustering of absorbing gas (Scannameco et al. 2006) 
or (c) mere chance coincidence of randomly distributed absorbers 
jRauch etal]|l997h . The core of the Lyor absorption in all three 
systems where (5v(Ovi) > 200 km s"' is highly saturated whereas 
the high velocity components show weak Lyo- absorption and rel- 
atively strong metal absorption. This could mean that the gas in 
these high velocity components is highly ionized and possibly of 
high metallicity. On the other hand, if the three systems where 
^v(Civ) > 200 km s"' follow the correlation seen in Fig. [H] we 
expect them to have (5i'(Ovi) > 200 km s"'. However, 5v(0vi) 
could not be measured due to either low S/N (in two cases) or 
blending with strong Lya lines. The system at Zabs= 2.2750 towards 
HE 2347-4342 is showing unsaturated well connected wide spread 
Lytt absorption with signature of high ionization. For the other two 
systems strong C iv is seen with heavily saturated Lya. The system 
at Zabs= 2.8265 towards HE 0940-1050 which shows largest veloc- 
ity spread {6v{C iv) ~ 375 km s"') for C iv in our sample is possibly 
associated with a Lyman break galaxy (LBG) at an impact parame- 
ter of 150 h"' kpc ( Crighton et al. 201 1). Hence it is extremely im- 
portant to have a detailed spectroscopic survey of galaxies around 
the redshifts of such large Sv systems to understand the possible 
origin of O vi absorbers. 



5.2 Correlation of 6v with other parameters 

In this section we explore any possible correlation between Sv and 
other observable parameters (see Fig. I13t . The strongest correla- 



tion is seen between 6v(0\i) and A'(Ovi). The Spearman rank 
correlation coefficient is 0.72 and the correlation is confirmed at 
5.3cr level. While systems with low 6v are seen over a wide range 
of A'(Ovi), the systems with <5v > 100 km s"' are seen when 
log ^(O vi) > 14. Thus the above mentioned correlation is due to 
the lower envelop one can see in the figure. As Sv is mainly re- 
lated to the number of components, this lower envelop can not be 
attributed to a detection bias. This is confirmed by the presence of a 
significant correlation (i.e. p, = 0.64 with a 4.0cr significance) even 
when we restrict our analysis to systems with log N{Owi) > 13.7. 

Interestingly, a similar correlation is also seen between 
N{Ci\) and <5v(Civ) (i.e. = 0.57 with a 5.8(T significance). The 
correlation is significant even when we restrict our analysis to sys- 
tems with log A'(Civ) > 12.6 (i.e. p, = 0.42 with a 3.7o" signifi- 
cance). These r esults are consistent with what is seen in Fig. 9 of 
'Song ailal JioO^ . who found C iv systems with peak optical depth 
greater than 0.4 to show large r velocity extent compared to systems 
with lower peak optical depth. lSongailal j2006h suggested that some 
of the wider C iv systems found among those with high peak optical 
depth could be associated with galaxy outflows. 

We also find a 3cr correlation between Sv{0 vi) and A'(H i) and 
a 4.5cr correlation between (5v(Civ) and A'(Hi). It is clear from 
Fig. 1 1 31 that for both Ovi and Civ, Sv > 100 km s"' is seen 
mainly in systems with log A'(Hi) > 15. Such high column den- 
sitie s are generally associated with high over density regions (see 
e.g.. lSchav3l200lh . In the following section, we show that the sys- 
tems with log Af(Hi) > 15.0 have associated low ion absorption. 
Low ions usually trace higher density regions if the gas is pho- 
toion ized by the meta-galactic UV radiation (see e.g.. Ranc h et aU 
Il997h . All these suggest that the O vi or C rv absorbing gas with 
high velocity spread is probably related to overdense regions. We 
do not find any strong correlation between Sv{0 vi) or Sv{C iv) with 
either redshift, ^(O vi)/A'(C iv) or ^(O vi)/A'(H i). 
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Figure 12. Absorption profiles of systems with Sv(0\i) > 200 km s"' (top panel) and Sv(Civ) > 200 km s"' (bottom panel) together with best fitted Voigt 
profile. The high velocity components identified with solid horizontal bars are showing high ionization. The absorption redshift that defines the zero velocity 
and the name of the background QSO are also indicated. 
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Figure 13. Relationship between Sv and other measurable quantities. The results of Spearman test performed for the systems using only the detections of Sv 
(filled circles) are mentioned in each panel. 



5.3 Velocity shift (Av) distribution 

In this section, we discuss the velocity olfset between the O vi and 
C IV absorption originating from the same system. Here, we do not 
attempt to estimate the offset between individual components as 
our main aim is to find the optical depth weighted phase separation 
between C iv and O vi and how this is related to other measurable 
quantities. In most of the cases, the O vi absorption is wider than the 
Civ one (see section [STt and has a larger number of components. 
We use only systems where both Civ and Ovi profiles are well 



defined at least in one of the doublets and calculate the optical depth 
weighted redshifts (i.e. zq yi and zq ly)- The velocity shift between 
these two redshifts is then, 

|Av(Ovi - Civ)| = [|(zo VI - 2Civ)l/(l + zCiv)] x (6) 

where c is the speed of light in km s"' . The values of |Av(0 vi-C iv)| 
measurements are summarized in Table[T] In the left most panel of 
Fig- El we plot the distribution of |Av(0 vi - C iv)|. The values of 
the velocity shift are found to be in the range < |Av(0 vi - C iv)| < 
48 km s"' with a median value of 8 km s"' . We find that only ~ 9% 
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Figure 14. Left : The distribution of the velocity shift between the optical depth weighted redshifts of C iv and O vi. Middle : The velocity shift, |Ai'(0 vi - C iv)| 
versus Sv{C iv) (bottom) and dv(0 vi) (top). Right : The velocity shift, |Av(0 vi - C iv)| versus ^(O vi) (bottom) and ^(H i) (top). The dashed lines in the middle 
and right panels are for illustrative purpose to emphasize the possible presence of a lower envelop. The (red) stars in both the panels represent systems with 
low ions. The results of Spearman rank congelation analysis for all the data points are also summarized in each panel. 



of the systems show |Av(Ovi - Crv)| > 20 km s"' whereas ~ 33% 
of the systems have |Av(Ovi - Civ)| < 5 km s"'. It is to be re- 
membered that the systems with large variation in A'(0 yi)/N(C iv) 
ratio from one component to another tend to have large values of 
|Av(0 VI - C rv)|. Hence |Ai'(0 vi - C iv)| can be seen as a measure 
of the ionization inhomogeneity among components. Our analysis 
suggests that Ovi systems with large ionization inhomogeneities 
across the profile are rare. In addition, the median value of the ve- 
locity shift distribution for systems with low ions (~ 10.4 km s"') 
is higher than that for systems without low ions (~ 6.8 km s"'). 
Hence, the ionization inhomogeneity seems to be more important 
in the systems where low ions are detected. But due to the small 
number of data points the KS test do not show any significant dif- 
ference between the shift distributions for systems with and without 
low ions. 

We find a mild (at ~ 2.5cr significance level) correlation be- 
tween either Sv(0 vi) or A''(0 vi) and |Av(0 vi - C iv)| (see Fig. 1141). 
These correlations are dominated by the fact that systems with 
|Ai'(0 VI - Civ)| > 20 km s"' are predominantly coming from sys- 
tems with 5v(0vi) > 100 km s"' and log A'(0 vi)(cm"-) > 14. 
A similar trend is also seen for (5v(Civ) and A'(Hi) possibly 
due to the presence of a lower envelop. However, we find none 
of the other parameters (i.e., z, N(Ovi)/N(lii), A'(Civ)/A'(Hi) 
and N(Ovi)/N(Crv)) showing any significant correlation with 
|Av(Ovi-Civ)|. 



6 ANALYSIS BASED ON TOTAL COLUMN DENSITIES 

Here we study the distributions of measured column densities of 
different species, their ratios and dependencies between them. Up- 
per limits are not considered for the analysis. 



6.1 Redshift evolution of column density ratios 

In Fig. [15] we plot the redshift evolution of various column density 
ratios. In all the panels (red) stars represent systems with detectable 
low ions and (black) open circles are for systems without low ions. 
The Spearman rank correlation analysis performed between 



various column density ratios and redshift do not reveal any sig- 
nificant correlation. We also find that the presence of low ions 
does not influence this result. One of the interesting features of 
the left most panel in Fig. [15] is the sc arcity of data points with 
log 7V(Ovi)/jV(Hi) > -0.5 for z > 2.5. [Bergeron & Herbert-For3 
(2005) identified such systems as a sepa rate population of O vi ab- 
sorbers with high metallicity (see also ISchaveetal.ll2007h . Such 
high value s of N(Owi)/N(Hi) ratio are also s een in proximate O vi 
absorbers jFox et alJlBoH : iTripp et alj |2()08l) . Therefore, these ab- 
sorbers may either trace high metallicity gas or regions ionized by 
AGN like sources. Hence, confirming the redshift evolution of these 
absorbers will be very important. It is to be noted that similar trend 
(i.e., lack of data points with high ^(C iv) to A'(H i) ratio for z > 2.5) 
is apparent from the right most panel of Fig.[T5] The median values 
of various column density ratios are also shown in Fig.[T5]with hor- 
izontal lines. The difference of median values for systems with and 
without low ions is maximum for ^(C iv)/N{0 vi) ratio (factor ~ 4). 
A two sided KS test shows that their distributions are different with 
very high significance (D = 0.55 & Prob. = 2.6x10"'). 

6.2 Relationship between column densities of various species 
and viability of photoionization model 

In this section, we investigate the total column densities of vari- 
ous species and any correlation between them. We compare these 
with simpl e photoionization m odel predictions using "CLOUDY 
v(07.02)" i lFerlandetaI.1 11998'). We would like to point out here 
that such a single phase photoionization model may oversimplify 
the problem in view of the possible multiphase structure of the 
absorbing gas, but such models are often used as a guideline to 
draw broad conclusions on the nature of the Ovi absorbers (e.g. 
see.isimcoe et al . 2002; Carswell et al. 2002; Bergeron et al. 200j; 
Bergeron & Herbert-Fort,, 2005. ; .Schave et aLi2007i ; .Muzahid et al.l 
20Ilh . 

The model assumes the absorbing gas to be an optically thin 
plane parallel slab illuminated by the meta-galactic UV back- 
ground radiation contributed by QSOs and Lyman break galax- 
ies at the median redshift (i.e. z = 2.3) of our survey. We use 
"HM05" background radiation available in CLOUDY v(07.02) 
based on the UV spectrum calculated using the method described in 
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Figure 15. Redshift evolution of O vi to H i (Left) Civ to O vi (Middle) and Civ to H i (Right) column density ratio. In all the panels the (red) stars are the 
systems with low ions and (black) open circles are the systems without any low ions. The vertical dotted line shows the median redshift of our sample. The 
median values of column density ratio for the full sample (solid horizontal line), systems with low ions (short dashed line) and systems without low ions (long 
dashed line) are also indicated in each panel. 



lHaardt&Madaiil ( ll996h . For simplicity, we assumed solar relative 
abundances and run our model with two different values of metal- 
licities (i.e. Z = 0.1 and 0.01 Zq). In this model, we incorporate the 
relation between t he particle den sity and the observed H i column 
density as given in lSchayj ( 1200 ll) . i.e., 



A'(Hi) 



2.7 X 10"cm-2(l + df'-T^"-^^ 



1+z 



QJi- 



0.02 i 1 0.16 



1/2 



X K, 



(V) 



where, 6 denotes over-density, fg is the gas mass fraction nor- 
malised to it's universal value (i.e. Q.h/ii„). All other symbols have 
their usual meaning. We have introduced a fudge factor called k 
and run this model for t hree different v alues of k (0.5, 1 and 2). For 
the original equation o f ISchavd JioOlh . K = I. Hence, for a given 
WCHi), K > I implies lower density ( higher ionizati on parame- 
ter) compared to what is predicted by ISchavd ( 1200 ih . The tem- 
perature is cal culated from the standard T - 6 relation given by 
iHui & Gnedinl (1997) i.e. T = To cP'"' with To = 2 x lO"* K and 
7 = 1.1. This temperature is very close to the photoionization tem- 
perature for the density range of interest here. 



6.2.1 N(Ovi)vs N(Hi) 

In the panel-(A) of Fig. [T6] we plot A'(0 vi) against A'(Hi) for indi- 
vidual O VI systems. Note that ^(O vi) is varying only ~ 2 dex over 
a ~ 5 dex variation in A'(Hi). The Spearman rank correlation anal- 
ysis shows only a mild correlation (at ~ 2.5cr level) between A'(Hi) 
and A'(0 vi) when we consider the full sample. The subsamples of 
systems with and without low ions do not show any trend between 
A'(H i) and ^(O vi) individually. 

It is apparent from the figure, that apart from two systems, low 
ions are detected in systems with higher Hi column density (i.e. 
log N(Hi) > 14.6). The median values of log A'(Hi) are 15.58 and 
14.82 for systems with and without low ions respectively. The KS 
test shows that the distributions of A'(Hi) in systems with and with- 
out low ions are different with very high significance (D = 0.52 and 
Prob. = 2.1x10"''). It is interesting to note that 36 out of 46 systems 
where we detect low ions show log A'(H i) > 15.0. This could mean 



that ^(H i) is predominantly coming from the low ionization phase 
of the absorbing gas whenever low ions are detected. This is why 
we do not use systems with low ions (even when metals are aligned 
with H i) when we measure thermal and non-thermal contribution 
to the /^-parameter in section l4~2l Although there is no statistically 
significant trend between A'(0 vi) and ^(H i), the O vi column den- 
sity seems to be systematically higher when low ions are present. 
The median value of log A'(0 vi) for systems with low ions is found 
to be 0.4 dex higher compared to that of the systems without low 
ions. A two sided KS test shows that the ^(O vi) distribution be- 
tween systems with and without low ions are indeed different (D = 
0.44 and Prob. = 7.3 x lO"*). 

It is evident from the model curves that the systems with 
log A'(Hi) > 14.5 can be roughly reproduced by our simple model 
for metallicity ranging from 0.01 to 0. 1 Zq. The systems in the top- 
right corner of the figure can also be reproduced with higher k val- 
ues (i.e. > 5) and with low metallicity (i.e. ~ O.OIZq). It is to be 
remembered that associating entire Hi column density measured 
in a system to the Ovi bearing phase will essentially underesti- 
mate the metallicity. If, in case of systems with low ions, most of 
the H I comes from the lower ionization phase then the (red) stars 
will move towards the left in this plot requiring higher metallic- 
ities. Most interestingly, the systems with log A'(Hi) < 14.5 and 
log A'(Ovi) > 14.0 can not be reproduced by our simple model 
with metallicity Z = O.OIZq (thin curves) irrespective of the k val- 
ues used. These systems require metallicity > 0. IZg provided the 
ionization is dominated by the meta-galactic UV background. 

In the low redshift studies, the A?(H i)/Af(0 vi) vs A'(Hi) plot 
is used as an indicator of multiphase medium. For collisional ion- 
ization, the Ovi ionization fraction peaks at ~ io^-''±"-^ K (see 
[Sutiierland & DoDitjl993ljGnat & Sternberd2007i ) and hence O vi 
is the most promising species in the UV-optical regime to probe the 
hard-to-detect hot gas phase in the intergalactic medium, namely, 
the WHIM. On the other hand most of the Lyo- lines are believed to 
trace relatively cool (T~ few x 10"* K) photoionized gas, namely, the 
warm ionized medium. To investigate the relative importance of the 
warm photoionized and hot collisionally ionized gas, it is custom- 
ary to plot the "multiphase ratio" i.e. A^(H i)/N (0 vi) against ^(H i) 
(see lShull et al.l2003| - |Danforth & Shulllioosl) . Panel-(B) of Fig.[ll] 




Figure 16. (A) O vi column density against H i column density. In all panels the (red) stars and the (black) circles represent systems with and without low ions 
respectively. The dot-dashed, solid, and long-dashed curves indicate the results of the photoionization model for k = 0.5, 1.0, and 2.0 respectively. The thin and 
thick curves correspond to metallicities O.OIZ© and 0. IZq. (B) Ovi to Hi colurnn density ratios in Ovi systems as a function of A'(Hi). The long-dashed line 
shows the power law fit to our data. The dotted line shows the fit obtained bv lDanforth & ShuU j2005i) for their low-z O vi sample. (C) C iv column density 
against H i column density. The shaded ones are the systems where only upper limits on N{0 vi) can be estimated. (D) A'(0 vi) IN{C iv) column density ratio 
against A'(Hi). Since the N(Ovi) to A'(Civ) column density ratio is independent of metallicity we show the curves for only Z = O.OIZq . 



is a very similar plot and it reveals a strong anti-correlation between 
A'(0 vi)/A'(Hi) and A'(Hi). This is not very surprising because it 
is the manifestation of ^(O vi) being not strongly correlated with 
A'(Hi). Note that the best fit straight line to our data (dashed line) 
gives a slope and interc ept very similar to what is obtained at low 
redshift (dotted line) bv' Panforth & Shulll feOOSi). We find that the 
systems with low ions show slightly higher A'(0 vi)/A'(Hi) value 
for a given A'(Hi). However, the slope of the best fitted straight 
lines for systems with and without low ions are very similar. 



6.2.2 N(Civ)vs N(Hi) 

Unlike A'(Ovi), N{Ciw) seems to be correlated with A'(Hi) (see 
panel-(C) of Fig. |16t. A Spearman rank correlation analysis sug- 
gests a correlation with rank coefficient ps = 0.41 and a 4.2cr sig- 
nificance level. We find that this correlation is mainly dominated by 
the systems with low ions. In fact no trend is seen between A'(C iv) 
and Af(Hi) for systems without low ions. The shaded points in this 
panel represent the systems where we have upper limits on Ovi 
column density as listed in Tabled and they are predominantly ly- 
ing in the lower half plane. The median value of log A^(C iv) (i.e. 
12.64) in these systems is 0.5 dex less than that for the rest of the 
systems. A two sided KS test indeed shows that the distribution of 
these points are different from the rest of the data points with ~ 98% 
probability. The non-detection of O vi in most of these systems is 
consistent with the average A'(0 \i)/N(C iv) ratio in our sample. 

The trend of increasing A'(Civ) with A'(Hi) is well repro- 
duced by our simple model. Unlike in the case of ^(O vi), the dif- 
ference in the predicted curves at a given A'(Hi) and metallicity 
is very small for the range of k used in our models (in particu- 
lar for 14 < log Af(Hi) < 16). Again here, most of the systems 
with log A'(H i) > 14.5 can be well accommodated with the model 
curves for metallicity ranging from 0.001 to 0.1 Zq. On the other 
hand systems with log A'(Hi) < 14.5 and log A'(Civ) > 13 require 
metallicity Z > 0. IZq irrespective of k. 



6.2.3 N(Ovi )/N(C IV ) vs N(Hi) 

In panel-(D) of Fig. [161 the ^(O vi}/N(C iv) ratio is plotted against 
the Hi column density. A strong anti-correlation (~ 4.3cr level) 



is seen between A'(Ovi)/A'(Civ) and A'(Hi). Since A'(Ovi) shows 
a mild correlation whereas A'(C iv) shows strong correlation with 
^(H i), such anti-correlation is expected. It is to be noted that the 
anti-correlation seen in the full sample is again dominated by the 
systems with low ions. No trend is suggested by Spearman rank 
correlation analysis for the systems without low ions. As the O vi 
to C IV column density ratio mainly depends on the ionization pa- 
rameter in this single phase model both the thick and thin curves 
fall on top of each other. Note that the observed anti-correlation is 
well reproduced by our model apart from a few systems with very 
low A'(Hi). 



7 SUMMARY & CONCLUSIONS 

We presented a detailed analysis of 84 O vi and 105 C iv systems 
at z ~ 2.3 detected in high resolution (R ~ 45,000) spectra of 18 
bright QSOs observed with VLT/UVES. Here we summarize the 
main results of our survey. 

• Multiphase nature of the gas : Consistent with the previous 
studies, we show that Ovi components have systematically wider 
Doppler parameters (b) compared to those of Civ components. We 
also show that the line spread (^v) of C iv and O vi are strongly cor- 
related. Therefore we conclude that the metal absorbing regions in 
the IGM consist of multiphase gas correlated over large velocities. 

We do not find any trend between A'(0 vi) and A'(H i), over five 
orders of magnitude spread in A'(H i) but there is a 3cr level corre- 
lation between line spread of O vi and A''(H i) which possibly sug- 
gests that the H i and O vi oc cur in different phases of a correlated 
structure. Indeed. Ifo^ ilO 1 ih recently argued that the constancy of 
^(O vi) over a range of ^(H i) can be reconciled if O vi absorp- 
tion originate from conductive, turbulent or shocked boundary lay- 
ers between warm and hot plasma. Such models are considered to 
explain the multiphase structure seen in different components of 
our galaxy (Savage & Sembach 1994; Savage et al. 2003; Spitzei 
199d ; ISembach et aTl2003l ; ICollins et all2004 12005| ; lUehneretal 



201 11) . Existing simulations of metals in the IGM roughly predict 
such a multiphase correlated structure where O vi absorption origi- 
nates from low density and extended region whereas C iv originates 
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from regio ns of slight ly higher density (e.g., iRauch et aljll997l: 
Kawata & R anch 2 007l:lFangano et alj2007l ; IOipp enheim er & Davd 
2009l : ICen & Chisar1l201lh . 



• Thermal state of the gas and non-thermal velocities : The ob- 
servations presented here are inconsistent with most of the gas asso- 
ciated with the O vi absorption (at z ~ 2.3) being in coUisional ion- 
ization equilibrium (i.e. with T ~ (2-3)x lO"" K). We draw this con- 
clusion mainly from the b(0 vi) distribution. Using a subsample of 
well aligned O vi components we find the median gas temperature 
is ~ 3 X 10"* K with none of the components having T > 2 x 10^ K. 
However, 42% of these well aligned components show 4.6 < log T 
< 5.0, which is warmer than the temperature expected in pure pho- 
toionization equilibrium. In case of rapidly cooling over-ionized 
gas such temperatures are expected provided the g as abundance is 
close to the solar value (see lGnat & Stemberd2007f) . The estimated 
non-thermal contribution to the fe-parameter using the fe(Ovi)- 
i>(Hi) pairs are in the range 3.6 < b„,(km s"') < 21.2 with a median 
value of 8.2 km s"'. This is consist ent with the previous mea sure- 
ments of intergalactic turbulence bv lRauch et al.l i [l996ll200lh . 

• Column density distribution : The Ovi column density dis- 
tribution is well fitted by a power-law with indices fj = 1.9 ± 
0.1 and 2.4 ± 0.2 for systems and components respectively for 
log A'(Ovi) > 13.7. We find the distribution is flatter in the case 
of Crv with respective values of 1.6 ± 0.1 and 1.9 ± 0.1 down 
to log N{Ci\) = 12.6. For both Civ and Ovi the /3 values mea- 
sured for the components are higher compared to that measured 
for systems. This is a natural consequence of the fact that the 
systems with higher column density having systematically higher 
number of components. While we could not make precise com- 
parisons of these results with the predictions of existing simula- 
tions, it is interesting to note that some of the simulations (with or 
wi thout feedback) pr o duce similar trends (se e for example Fig. 10 
of lRauch et Zlll997h . ICen & ChisMil ( 1201 ih have found a steeper 
slope for the N(0 vi) distribution compared to that of C iv. They at- 
tributed this to the existence of transient structures stemming from 
the shock heated regions in the neighborhood of galaxies. 

• Gas kinematics : A strong correlation (5.3cr level) is seen be- 
tween the velocity spreads (6v) of O vi and C iv. However, 6v(C iv) 
is systematically lower than (5v(0vi). We find that both 6v{Ovi) 
and Sv(C iv) are strongly correlated (> 5<t level) with their respec- 
tive column densities and slightly less (3 and 4.5(T respectively) 
correlated with A'(H i). We note that O vi and/or C iv systems with 
large velocity spread also show associated low ion absorption lines. 
As the low ions (as well as strong H i absorption) are expected to 
originate from high density regions, we conclude that systems with 
large velocity spread are probably associated with regions of high 
density. We measured the velocity offset, |Av(0 vi - C iv)|, between 
optical depth weighted redshifts of C iv and O vi absorption, which 
is found to be in the range < |Av(Ovi - Civ)| < 48 km s"' with 
a median value of 8 km s"'. The systems with low ions seem to 
have higher velocity shift which possibly indicates higher ioniza- 
tion inhomogeneity in these absorbers. We do not find any strong 
correlation between |Av(0 vi - C iv)| and other observable parame- 
ters. 

• Total column densities and their ratios : The total column 
densities of different species (i.e. Hi, Civ, and Ovi) seem to be af- 
fected by the presence of low ions. The median values of A'(H i), 
A^(Ovi) and A'(Civ) are found to be higher when low ions are 
present. A two sided KS test suggests that the column density 



distributions for systems with and without low ions are signif- 
icantly different. Almost ~ 80% of the systems with low ions 
show log A'(Hi) > 15.0 indicating that considerable Hi absorp- 
tion may be originating from the low ionization phase. We find 
a strong correlation (~ 4.3cr) between A'(Civ) and A'(Hi) which 
is dominated by the systems with low ions. We do not find any 
clear evidence for the column density ratios (i.e. ^(O vi)/A'(Hi), 
A'(Civ)/A'(Hi) and ^(O vi)/A'(Civ)) to evolve with redshift over 
the range 1.9 < z < 3.1. We find only a tentative evidence for num- 
ber of systems with high log A'(0 vi)/Af(Hi) (e.g., > -0.5 dex) to 
decrease with increasing redshift. A similar trend is also present for 
^(C iv)/A'(Hi) ratio. 

• Comparison between low and high-z O vi absorbers : The 

/7-parameter distribution of O vi components of our sa mple is sig- 
nificantly different from that of the low-z sample of iTripp et al.l 
(2008). The medi an value o f fc(Ovi) at low-z is twice as high as 
at high-z (see also lFoxl|201 ih . In addition, the median fc„t value for 
the well aligned components in our sample is a factor ~ 2 less than 
what is found bv lTripp et alj ( 1200 8l) . Interestingly the median value 
of the temperature measured in these well aligned components is 
found to be the same (i.e. T ~ 3 x 10* K) in both high and low-z 
sample. All these are consistent with the non-ther mal contribution 
to b(0 vi) being higher at low-z. In the models of lEvoli & Ferraral 
( [201 l|) where the turbulence induced by supernova driven winds 
is considered, the volume-weighted i-parameter remains roughly 
constant between z = 2 to z = [see their Fig. 6]. Such a case 
is not supported by our observations. We speculate that the excess 
of turbulence at low-z may originate from shocks due to structure 
formation. 

At low redshift, a few thermally broadened Lyo- absorbers 
(BLAs) are seen with signature of high temperature ( T ~ 10^ K) 
gas ( Sembach et al. 2004; Richter et al. 2004, 2006; D anforth et al] 
1 20 id ; ISavage et alj|2011h . We do not detect such systems in our 
high-z sample. This may be related to difficulties in detecting broad 
and shallow absorption features in the dense Lya forest. 

At z ~ 2.3 we find the cosmic density of O vi absorbers, Q.q yj 
= (I.0±0.2)xl0"' for log A^(0 vi) > 13.7. This should be treated as 
a lower limit as (a) we have not applied the correction factor to red- 
shift path length due to Lyo- line blanketing in the forest, (b) the 
CDDF of O VI is steep and Q.q yj may increase when contributions 
of numerous low column density systems are included. We calcu- 
late the lower limits on the baryonic content of the O vi absorbers 
assuming (i) ionization fraction of O vi, /q yj = 0.2, (ii) the average 



sumed in the low-z studies of vi absorbers (see, Tripp et al. 


200C 


Savage et al. 


2002 


; Sembach et alj 2004; Danforth & Shulll 


2005 


Lehner et al 


[2006 


; [Danforth & Shullll2008l). Most of these papers 



have shown that low redshift O vi absorbers harbor roughly 5-10% 
of the baryons in the nearby universe. We find this contribution to 
be 2.8% at z ~ 2.3. The correcti on to the redshift pa th length due to 
line blanketing (as suggested bv [Simcoe et alj2002h will increase it 
up to ~ 7%. Therefore within allowed uncertainties, flj^j^' at z ~ 
is consistent with what we found at z ~ 2.3. If Ovi predominantly 
traces a hot phase of the IGM (i.e. T > 10^ K) at every epoch 
then most numerical simulations ( Cen & Ostrikeil 19991; Dave et all 
[200 ll ; [Fang & Brvan[|200ll ; IChen et alj|2003l ; ICen & Chisarill20Ilh 
suggest a strong increase in f^f^J' with decreasing redshift. These 
simulations also suggest that the fraction of baryons at T < 10' K 
decreases with decreasing z. Therefore the near constancy of fi^jj' 
with redshift probably means that the Ovi absorbers at high and 



20 S. Muzahid et al. 



low redshift may not originate from regions with similar physical 
conditions. 
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